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CO Abstract 

<^ Leptogenesis is a class of scenarios in which the cosmic baryon asymmetry originates from an initial 
^ lepton asymmetry generated in the decays of heavy sterile neutrinos in the early Universe. We ex- 
• • plain why leptogenesis is an appealing mechanism for baryogenesis. We review its motivations, the 
basic ingredients, and describe subclasses of effects, like those of lepton flavours, spectator processes, 
scatterings, finite temperature corrections, the role of the heavier sterile neutrinos and quantum cor- 
^ rections. We then address leptogenesis in supersymmetric scenarios, as well as some other popular 
variations of the basic leptogenesis framework. 
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1.1 The Baryon Asymmetry of the Universe 
1.1.1 Observations 

Up to date no traces of cosmological antimatter have been observed. The presence of a small amount 
of antiprotons and positrons in cosmic rays can be consistently explained by their secondary origin 
in cosmic particles collisions or in highly energetic astrophysical processes, but no antinuclei, even as 
light as anti-deuterium or as tightly bounded as anti-a particles, has ever been detected. 

The absence of annihilation radiation —)•... vr'^ —)•... 27 excludes significant matter-antimatter 
admixtures in objects up to the size of galactic clusters ~ 20Mpc f^. Observational limits on anoma- 
lous contributions to the cosmic diffuse 7-ray background and the absence of distortions in the cosmic 
microwave background (CMB) implies that little antimatter is to be found within ~ 1 Gpc and that 
within our horizon an equal amount of matter and antimatter can be excluded [2]. Of course, at larger 
super- horizon scales the vanishing of the average asymmetry cannot be ruled out, and this would 
indeed be the case if the fundamental Lagrangian is C and CP symmetric and charge invariance is 
broken spontaneously [3]. 

Quantitatively, the value of baryon asymmetry of the Universe is inferred from observations in two 
independent ways. The first way is by confronting the abundances of the light elements, D, ^He, ^He, 
and ^Li, with the predictions of Big Bang nucleosynthesis (BBN) [H O [6l [3 [8l |9] . The crucial time for 
primordial nucleosynthesis is when the thermal bath temperature falls below T ^ 1 MeV. With the 
assumption of only three light neutrinos, these predictions depend on a single parameter, that is the 
difference between the number of baryons and anti-baryons normalized to the number of photons: 



where the subscript means "at present time" . By using only the abundance of deuterium, that is 
particularly sensitive to r/, Ref. 0] quotes: 

10^% = 5.7 ±0.6 (95%c.l.). (1.2) 

In this same range there is also an acceptable agreement among the various abundances, once theo- 
retical uncertainties as well as statistical and systematic errors are accounted for [6J. 

The second way is from measurements of the CMB anisotropics (for pedagogical reviews, see 
Refs. |10| nrij). The crucial time for CMB is that of recombination, when the temperature dropped 
down to T ;^ 1 eV and neutral hydrogen can be formed. CMB observations measure the relative 
baryon contribution to the energy density of the Universe multiplied by the square of the (reduced) 
Hubble constant h = i7o/(100km sec"^ Mpc"^): 

^sh'^^h^^, (1.3) 

Pcrit 

that is related to r] through lO^"?? = 274 Q^b h'^- The physical effect of the baryons at the onset of matter 
domination, which occurs quite close to the recombination epoch, is to provide extra gravity which 
enhances the compression into potential wells. The consequence is enhancement of the compressional 
phases which translates into enhancement of the odd peaks in the spectrum. Thus, a measurement of 
the odd/even peak disparity constrains the baryon energy density. A fit to the most recent observations 
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(WMAP7 data only, assuming a ACDM model with a scale-free power spectrum for the primordial 
density fluctuations) gives at 68% c.l. |T2] 



102QB/i' = 2.258l°:°f^. (1.4) 

There is a third way to express the baryon asymmetry of the Universe, that is by normalizing the 
baryon asymmetry to the entropy density s = 5*(27r^/45)r^, where is the number of degrees of 
freedom in the plasma, and T is the temperature: 



(1.5) 

The relation with the previous definitions is given by the conversion factor sq/u^q = 7.04. Ias is a 
convenient quantity in theoretical studies of the generation of the baryon asymmetry from very early 
times, because it is conserved throughout the thermal evolution of the Universe. In terms of Iab the 



BBN results eq. (1.2) and the CMB measurement eq. (1.4) (at 95%c./.) read: 



Y^B^ = (8-10 ± 0.85) X 10"^\ Y^^'^ = (8.79 ± 0.44) x 10"^^ (1.6) 

The impressive consistency between the determinations of the baryon density of the Universe from 
BBN and CMB that, besides being completely independent, also refer to epochs with a six orders of 
magnitude difference in temperature, provides a striking confirmation of the hot Big Bang cosmology. 



1.1.2 Theory 

From the theoretical point of view, the question is where the Universe baryon asymmetry comes from. 
The infiationary cosmological model excludes the possibility of a fine tuned initial condition, and since 
we do not know any other way to construct a consistent cosmology without inflation, this is a strong 
veto. 

The alternative possibility is that the Universe baryon asymmetry is generated dynamically, a 
scenario that is known as baryogenesis. This requires that baryon number (B) is not conserved. More 
precisely, as Sakharov pointed out [l3], the ingredients required for baryogenesis are three: 

1. B violation is required to evolve from an initial state with Iab = to a state with 1a_b 7^ 0. 

2. C and CP violation: If either C or CP were conserved, then processes involving baryons would 
proceed at the same rate as the C- or CP-conjugate processes involving antibaryons, with the 
overall effect that no baryon asymmetry is generated. 

3. Out of equilibrium dynamics: Equilibrium distribution functions ricq are determined solely by 
the particle energy E, chemical potential /x, and by its mass which, because of the CPT theorem, 
is the same for particles and antiparticles. When charges (such as B) are not conserved, the 
corresponding chemical potentials vanish, and thus hb = J (^^'^eq = ng. 

Although these ingredients are all present in the Standard Model (SM), so far all attempts to reproduce 
quantitatively the observed baryon asymmetry have failed. 

1. In the SM B is violated by the triangle anomaly. Although at zero temperature B violating 
processes are too suppressed to have any observable effect ^4j , at high temperatures they occur 
with unsuppressed rates |15j . The flrst condition is then quantitatively realized in the early 
Universe. 
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2. SM weak interactions violate C maximally. However, the amount of CP violation from the 
Kobayashi-Maskawa complex phase [16], as quantified by means of the Jarlskog invariant |17j. is 
only of order 10"^*^, and this renders impossible generating Iab ~ 10"^*^ |18 | 1191 [20] . 

3. Departures from thermal equilibrium occur in the SM at the electroweak phase transition 
(EWPT) |2H I22| . However, the experimental lower bound on the Higgs mass implies that 
this transition is not sufficiently first order as required for successful baryogenesis [23j. 

This shows that baryogenesis requires new physics that extends the SM in at least two ways: It must 
introduce new sources of CP violation and it must either provide a departure from thermal equilib- 
rium in addition to the EWPT or modify the EWPT itself. In the past thirty years or so, several 
new physics mechanisms for baryogenesis have been put forth. Some among the most studied are 
GUT baryogenesis [Ml[2a[Ml[2Zl[28l[29l|30l[3ll|32l[33l Electroweak baryogenesis [2ll [Ml ES] , the 
Affleck-Dine mechanism [361 ET], Spontaneous Baryogenesis [3Sll39]. However, soon after the discov- 
ery of neutrino masses, because of its connections with the seesaw model jlQl HH HJl |33l 01] and its 
deep interrelations with neutrino physics in general, the mechanism of baryogenesis via Leptogene- 
sis acquired a continuously increasing popularity. Leptogenesis was first proposed by Fukugita and 
Yanagida in Ref. [l5] . Its simplest and theoretically best motivated realization is precisely within the 
seesaw mechanism. To implement the seesaw, new Majorana SU{2)l singlet neutrinos with a large 
mass scale M are added to the SM particle spectrum. The complex Yukawa couplings of these new 
particles provide new sources of CP violation, departure from thermal equilibrium can occur if their 
lifetime is not much shorter than the age of the Universe when T ~ M, and their Majorana masses 
imply that lepton number is not conserved. A lepton asymmetry can then be generated dynamically, 
and SM sphalerons will partially convert it into a baryon asymmetry |l6]. A particularly interesting 
possibility is "thermal leptogenesis" where the heavy Majorana neutrinos are produced by scatterings 
in the thermal bath starting from a vanishing initial abundance, so that their number density can be 
calculated solely in terms of the seesaw parameters and of the reheat temperature of the Universe. 



This review is organized as follows: in Section 1.2 the basis of leptogenesis are reviewed in the 
simple scenario of the one fiavour regime, while the role of fiavour effects is described in Section 1.3 
Theoretical improvements of the basic pictures, like spectator effects, scatterings and CP violation in 
scatterings, thermal corrections, the possible role of the heavier singlet neutrinos, and quantum effects 



are reviewed in Section 1.4 Leptogenesis in the supersymmetric seesaw is reviewed in Section 1.5 



while in Section 1.6 we mention possible leptogenesis realizations that go beyond the type-I seesaw. 



Finally, in Section 1.7 we draw the conclusions. 



1.2 Ni Leptogenesis in the Single Flavour Regime 



The aim of this section is to give a pedagogical introduction to leptogenesis |45j and establish the 
notations. We will consider the classic example of leptogenesis from the lightest right-handed (RH) 
neutrino A''i (the so-called A^i leptogenesis) in the type-I seesaw model [lOl 031 HH HI] in the single 



fiavour regime. First in Section 1.2.1 we introduce the type-I seesaw Lagrangian and the relevant 



parameters. In Section 1.2.2, we will review the CP violation in RH neutrino decays induced at 1-loop 



level. Then in Section 1.2.3, we will write down the classical Boltzmann equations taking into account 
of only decays and inverse decays of A'^i and give a simple but rather accurate analytical estimate of 
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the solution. In Section 1.2.4 we will relate the lepton asymmetry generated to the baryon asymmetry 



of the Universe. Finally in Section 1.2.5 we will discuss the lower bound on A^i mass and the upper 



bound on light neutrino mass scale from successful leptogenesis. 
1.2.1 Type-I seesaw, neutrino masses and leptogenesis 

With m (m > 2^^ singlet RH neutrinos Nj^. (i = 1, m), we can add the following Standard Model (SM) 
gauge invariant terms to the SM Lagrangian 

Ci = CsM + iW.0NR^ - (^^MiN^^NR^ + CabYaiW^r^H' + h.c}j , (1.7) 

where Mj are the Majorana masses of the RH neutrinos, ia = {^aLidJ,) with a = e, ^, r and H = 
{H^, H^) are respectively the left-handed (LH) lepton and Higgs SU{2)l doublets and e^b = —e^a with 
ei2 = 1. Without loss of generality, we have chosen the basis where the Majorana mass term is diagonal. 
The physical mass eigenstates of the RH neutrinos are the Majorana neutrinos Ni = Nr. + Nf^, . Since 
Ni are gauge singlets, the scale of Mi is naturally much larger than the electroweak (EW) scale 
Mi ^ (<I>) = V = 174 GeV. Hence after EW symmetry breaking, the light neutrino mass matrix is 
given by the famous seesaw relation jlQl HSj HU IB] 

~ -v^Y^Y^. (1.8) 



Assuming Y ~ 0(1) and rrii, ~ yAm^^^ ~ 0.05 eV, we have M ~ 10^^ GeV not far below the GUT 
scale. 

Besides giving a natural explanation of the light neutrino masses, there is another bonus: the three 



Sakharov's conditions |jl3j for leptogenesis are implicit in eq. (1.7) with the lepton number violation 
provided by Mi, the CP-violation from the complexity of Yia and the departure from thermal equilib- 
rium condition given by an additional requirement that Ni decay rate T]^^ is not very fast compared 
to the Hubble expansion rate of the Universe H{T) at temperature T = Mi with 

^N- = - H(T) = -\- , 1.9 

where Mpi = 1.22 x 10^^ GeV is the Planck mass, 5=1, (=106.75 for the SM excluding RH neutrinos) is 
the total number of relativistic degrees of freedom contributing to the energy density of the Universe. 
To quantify the departure from thermal equilibrium, we define the decay parameter as follows 

where irn is the effective neutrino mass defined as|47] 



(yty)..t,2 



with m=K = "fff^ Y ~ 1 X 10 ^ eV. The regimes where Ki <^ 1, Ki ^ 1 and Ki ^ 1 are respectively 
known as weak, intermediate, and strong washout regimes. 

^Neutrino oscillation data and leptogenesis both require m > 2. 
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1.2.2 CP asymmetry 

The CP asymmetry in the decays of RH neutrinos Ni can be defined as 

7 (iV, ^ iaH) --f[Ni^ I^H*) 

where 7(i — )■ /) is the thermally averaged decay rate defined a^ 



Ay 

IN, 



a 
Ni 



(1.12) 



-(27r)V^)(p,-p^)|^(z^/)|V 



where A{i 



f) is the decay amplitude. Ignoring all thermal efi'ects 



(1.13) 



Ea lA {Ni ^ laH) |2 + \Ao {Ni ^ £aH^ 



|9], eq. (1.12) simplifies to 
(1-14) 



where Ao{i — )• /) denotes the decay amplitude at zero temperature. Eq. ( 1.14 ) vanishes at tree level but 
is induced at 1-loop level through the interference between tree and 1-loop diagrams shown in Figure 
|1.1[ There are two types of contributions from the 1-loop diagrams: the self-energy or wave diagram 
(middle) [50] and the vertex diagram (right) [l5]. At leading order, we obtain the CP asymmetry |51] : 



— , , , y Im [(yty) ..y^.v* 




where the loop function is 



9{x) 



-Fl - (1 -h x) In 



1 + x 

X 



(1.15) 



(1.16) 



The first term in eq. (1.15) comes from L- violating wave and vertex diagrams, while the second term 
is from the L-conserving wave diagram. The terms of the form (Mf — M?)~^ in eq. (1.15) are from the 

Mj (resonant 

leptogenesis scenario, see Section 



1.6.1 



otherwise the CP asymmetry vanishes because the Yukawa couplings combination becomes real. 



In the one fiavour regime, we sum over the fiavour index a in eq. (1.15) and obtain 

1 

0+^ 



1 

8^ 




(1,17) 



where the second term in eq. (1.15) vanishes because the combination of the Yukawa couplings is real. 



^Here the Pauli-blocking and Bose-enhancement statistical factors have been ignored and we also assume Maxwell- 
Boltzmann distribution for the particle i i.e. /i = e"^''^^. See Refs. [481 149j for detailed studies of their effects. 

■^Notice that the resonant term becomes singular in the degenerate limit Mi = Mj. This singularity can be regulated 
by using for example an effective field-theoretical approach based on resummation [52] . 
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Figure 1.1: The CP asymmetry in type-I seesaw leptogenesis results from the interference between tree 
and 1-loop wave and vertex diagrams. For the 1-loop wave diagram, there is an additional contribution 
from L-conserving diagram to the CP asymmetry which vanishes when summing over lepton flavours. 

1.2.3 Classical Boltzmann equations 

We work in the one flavour regime and consider only the decays and inverse decays of A^^i . If leptogenesis 
occurs at T ^ 10^^ GeV, then the charged lepton Yukawa interactions are out of equilibrium, and this 
defines the one flavour regime. The assumption that only the dynamics of A'^i is relevant can be 
realized if for example the reheating temperature after inflation is Tr// <^ ^2,3 such that are 
not produced. In order to scale out the effect of the expansion of the Universe, we will introduce the 
abundances, i.e. the ratios of the particle densities Ui = J d^pfi to the entropy density s = ^^g^T^: 

^ ^. (1.18) 

s 

The evolution of the A^i density and the lepton asymmetry Ial = SYa^ = 2(y£ — Y^)|^can be described 
by the following classical Boltzmann equations (BE)[53j 

^ = -D,{Y^,-Y;,1), (1.19) 

^ = eiD,{Y^,-Y^l)-WiYAL, (1.20) 
where z = Mi/T and the decay and washout terms are respectively given by 

Di{z) = J^^=K.z^^, WAz) = -B,{z)^, (1.21) 

with K,n the n-th order modified Bessel function of second kind. Y^^ and Y^'^ readj^ 



From eq. (1.19) and eq. (1.20), the solution for Yal can be written down as follows 



>al(2;) = Y/s,L{zi)e - j dz ei{z)—j-j-e J^' ' (1.23) 

where Zj is some initial temperature when A^i leptogenesis begins, and we have assumed that any 
preexisting lepton asymmetry vanishes Y^^{zi) = 0. Notice that ignoring thermal effects, the CP 



asymmetry is independent of the temperature ei{z) = ei (c.f. eq. (1.17)). 



*The factor of 2 comes from the SU{2)l degrees of freedoms. 

^To write down a simple analytic expression for the equilibrium density of A^'i, we assume Maxwell-Boltzmann distri- 

I 



bution. However, foUwing [S3], the normalization factor Yf'' is obtained from a Fermi-Dirac distribution. 



8 



Neutrino Physics 



Leptogenesis in the Universe 



Weak washout regime 

In the weak washout regime {Ki <C 1), the initial condition on the A^i density lAr^(zj) is important 



If we assume thermal initial abundance of A''! i.e. YAr^(zi) = Y^'^(O), we can ignore the washout when 



A^i starts decaying at z ^ 1 and we have 



YLioc) -61 dz'^ = e.Y^liO). (1.24) 

On the other hand, if we have zero initial A'^i abundance i.e. Y/v^(zj) = 0, we have to consider the 
opposite sign contributions to lepton asymmetry from the inverse decays when A^i is being populated 
(Yat^ < Y/^^) and from the period when A''! starts decaying (Y/Vi > ^^Vi)- Taking this into account the 
term which survives the partial cancellations are given by [55j 

Yl^{oo)c.'^e,K!Y^l{0). (1.25) 

Strong washout regime 

In the strong washout regime {Ki ^ 1) any lepton asymmetry generated during the A'^i creation phase 
is efficiently washed out. Here we adopt the strong washout balance approximation [56] which states 
that in the strong washout regime, the lepton asymmetry at each instant takes the value that enforces 
a balance between the production and the destruction rates of the asymmetry. Equating the decay 



and washout terms in eq. (1.20), we have 



" -W{lf'^ " -Wiz)''^ = ^^^^^ ' ^'-''^ 

where in the second approximation, we assume Yn-^ ~ ^Ni- The approximation no longer holds when 
Yal freezes and this happens when the washout decouples at zj i.e. W{zf) < 1. Hence, the final 
lepton asymmetry is given bj|^ 

V-,,(oo) = -|_,,y;,._ZL,,,..«,o). (1,27) 

The freeze out temperature Zf depends mildly on Ki. For Ki = 10- 100 we have for example Zf = 7- 
10. We also see that independently of initial conditions, in the strong regime Yal{oo) goes as K^^ . 

1.2.4 Baryon asymmetry from EW sphaleron 

The final lepton asymmetry Yal{oo) can be conveniently parametrized as follows 

i^AL(oo) = eir?iy^^(0), (1.28) 



where rji is known as the efficiency factor. In the weak washout regime {Ki <C 1) from eq. (1.24) 



we have r/i = 1 (= j^Kf < 1) for thermal (zero) initial A'^i abundance. In the strong washout regime 



{Ki » 1), from eq. (1.27), we have ryi = gjy^ < 1- 



This differs from the efflciency in Ref. [55j by the factor which is due to the different normalization eq. ( 1.22 1 
^Compare this to a more precise analytical approximation in Ref. [55) . 
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If leptogenesis ends before EW sphaleron processes become active {T ^ lO^^ Qgy)^ the S - L 
asymmetry Ias-l is simply given by 

Ya^^, = -Yal. (1.29) 

At the later stage, the B — L asymmetry is partially transferee! to a S asymmetry by the EW sphaleron 
processes through the relation |57j 



28 

Yab{oo) = — 1ab-l(c») , 



(1.30) 



that holds if sphalerons decouple before EWPT. This relation will change if the EW sphaleron processes 
decouple after the EWPT \57\ |58] or if threshold effects for heavy particles like the top quark and 
Higgs are taken into account [581 EH] ■ 

1.2.5 Davidson-Ibarra bound 

Assuming a hierarchical spectrum of the RH neutrinos (Mi <^ M2, M3) and that the dominant lepton 
asymmetry is from the A^i decays, from eq. (1.17) the CP asymmetry from A'^i decays can be written 
as 



ei 



Ml 
M 



(1.31) 



Assuming three generations of RH neutrinos (n = 3) and using the Casas-Ibarra parametrization [60] 
for the Yukawa couplings 

1 



Y 



(1.32) 



miv 



diag ( Ml , M2 , M3 ) , Dm,, = diag(mi,^, m-jyj, Jrij^s) and R any complex orthogonal matrix 



where 

satisfying R = RR^ = 1, eq. (|1.31[) becomes 



ei 



3 MiY 



167rt;2 Y,m.^\Ru\' 

i 

Using the orthogonality condition = 1, we then obtain the Davidson-Ibarra (DI) bound [61 

i 

3 Ml 



(1.33) 



DI 



167r V' 



3 Ml Aml,^ 



(1.34) 



where m,y^ {''^^ui) is the heaviest (lightest) light neutrino mass. Applying the DI bound on eqs. (1.28)- 
(1.30), and requiring that Yab{oo) > Y§^^^ ~ 10~^°, we obtain 



Ml 



0.1 eV 



^3 



7?r''(A^i) ^ 10^ GeV, 



(1.35) 



where the 7/j""^(Mi) is the efficiency factor maximized with respect to Ki eq. (1.10) for a particular 
value of Ml. This allows us to make a plot of region which satisfies eq. (1.35) on the (Mi,m,yj) 
plane and hence obtain bounds on Mi and rriy-^. Many careful numerical studies have been carried 
out and it was found that successful leptogenesis with a hierarchical spectrum of the RH neutrinos 
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requires Mi ^ 10^ GeV 121 163] and m^^ ^ 0.1 eV [HI ESI ESI |66]. This bound impUes that the 
RH neutrinos must be produced at temperatures T ^ 10^ GeV which in turn imphes the reheating 
temperature after inflation has to be T^h ^ 10^ GeV in order to have sufficient RH neutrinos in the 
thermal bath. To conclude this section, let us note that the DI bound eq. (1.34) holds if and only if 
all the following conditions apply: 

(1) A^i dominates the contribution to leptogenesis. 

(2) The mass spectrum of RH neutrinos are hierarchical Mi <C M2, M3. 

(3) Leptogenesis occurs in the unflavoured regime T ^ 10^^ GeV. 

As we will see in the following sections, violation of one or more of the above conditions allows us to 
lower somewhat the scale of leptogenesis. 



1.3 Lepton Flavour Effects 

1.3.1 When are lepton flavour effects relevant? 

The first leptogenesis calculations were performed in the single lepton flavour regime. In short, this 
amounts to assuming that the leptons and antileptons which couple to the lightest RH neutrino 
Ni maintain their coherence as flavour superpositions throughout the leptogenesis era, that is £1 = 

Cqi^q and Ji = ^aC*]Ia- Note that at the tree-level the coefficients c and c'* are simply the 
Yukawa couplings: Cai = Yai and c'*i = Y*i. However it should be kept in mind that since CP is 
violated by loops, beyond the tree level approximation the antilepton state is not the CP conjugate 
of the £1, that is c'„^ 7^ Cai- 

The single flavour regime is realized only at very high temperatures (T ^ 10^^ GeV) when both 
ii and i'l remain coherent flavour superpositions, and thus are the correct states to describe the 
dynamics of leptogenesis. However, at lower temperatures scatterings induced by the charged lepton 
Yukawa couplings occur at a sufficiently fast pace to distinguish the different lepton flavours, ii and 
i'l decohere in their flavour components, and the dynamics of leptogenesis must then be described 
in terms of the flavour eigenstates ia- Of course, there is great interest to extend the validity of 
quantitative leptogenesis studies also at lower scale T ^ 10^^ GeV, and this requires accounting for 
flavour effects. The role of lepton flavour in leptogenesis was first discussed in Ref. [67 ], however 
the authors did not highlight in what the results were significantly different from the single flavour 
approximation. Therefore, until the importance of flavour effects was fully clarified in Refs. |68 |I691 [70] 
they had been included in leptogenesis studies only in a few cases [7T1 [721 E31 EH ES] . Nowadays lepton 
flavour effects have been investigated in full detail |76llZ7l|78l|79l[80l[8ll[8a[8a[8l[85llM^ 
and are a mandatory ingredient of any reliable analysis of leptogenesis. 

The speciflc temperature when leptogenesis becomes sensitive to lepton flavour dynamics can be 
estimated by requiring that the rates of processes Ta (a = e, /i, r) that are induced by the charged 
lepton Yukawa couplings ha become faster than the Universe expansion rate H(T). An approximate 
relation gives [90^ ,91j 

Ta{T) ^IQ-^lT , (1.36) 

which implies thal|^ 

r„(r) > H{T) when T ^ r„ , (1.37) 

*In supersymmetric case, since — rria/ {vu cos fi), we have T Tq (1 + tan^ /3). 
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where Te ~ 4 x 10^ GeV, T,, ~ 2 x 10^ GeV, and ~ 5 x 10^^ GeV. Notice that to fully distinguish 
the three flavours it is sufficient that the r and /i Yukawa reactions attain thermal equilibrium. It 
has been pointed out that besides being faster than the expansion of the Universe, the charged lepton 
Yukawa interactions should also be faster than the A'^i interactions [69\ [83l . In general whenever 
Tt-{Mi) > H{Mi) we also have Tr{Mi) > Tn-^(Mi). However, there exists parameter space where 
r^(Mi) > H{Mi) but Tr{Mi) < Tni{Mi). This scenario was studied in Ref. [83]. 



1.3.2 The effects on CP asymmetry and washout 

The CP violation in Ni decays can manifest itself in two ways [69j : 

(i) The leptons and antileptons are produced at different rates, 

where 7, = -/{Ni ^ £iH) and 7^ = -f{Ni ^ I,H*). 

(ii) The leptons and antileptons produced are not CP conjugate states, 

CP(£[) = / li, 

that is, due to loops effects they are slightly misaligned in flavour space. 



(1.38) 



(1.39) 



We can rewrite the CP asymmetry for Ni decays from eq. (1.12) as follows 



7i + 7i 



Pia ~l~ Pia , Pia 

Ci H 



Pu 



APi, 



(1.40) 



where terms of order 0{€i APia) and higher have been neglected. Pia is the projector from state ii 
into flavour state ia and APia = Pia — Pia- At tree level, clearly, Pia = Pia = Pj^a where the tree level 
flavour projector is given by 



(Y^Y)ii 



(1.41) 



Prom eq. ( 1.40 ), we can identify the two types of CP violation, the first term being of type (i) eq. (1.38 ) 
while the second being of type (ii) eq. (1.39). Since J2a Pia = J2a Pia = 1; when summing over flavour 
indices a, the second term vanishes J2a APia = 0. Note that the lepton- flavour- violating but L- 
conserving terms in the second line of eq. (1.15) is part of type (ii). In fact, they come from d = 6 
L-conserving operators which have nothing to do with the unique d = 5 L-violating operator (the 
Weinberg operator |92j) responsible for neutrino masses. However, in some cases they can still domi- 
nate the CP asymmetries but, as we will see in Section 1.3.4 lepton flavour equilibration effects [93] 



then impose important constraints on their overall effects. Note also that due to flavour misalignment, 
the CP asymmetry in a particular flavour direction €ia can be much larger and even of opposite sign 
from the total CP asymmetry ej . In fact the relevance of CP violation of type (ii) in the flavour regimes 



is what allows to evade the DI bound eq. (1.34). As regards the washout of the lepton asymmetry of 
flavour a, it is proportional to 

Wia (X Pia7^ + PiaJi ^ PiaWi, (1.42) 



which results in a reduction of washout by a factor of P^*^ < 1 compared to unflavoured case. As we 
will see next, the new CP- violating sources from flavour effects and the reduction in the washout could 
result in great enhancement of the final lepton asymmetry and, as was first pointed out in Ref. [69], 
leptogenesis with a vanishing total CP asymmetry ej = also becomes possible. 
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1.3.3 Classical flavoured Boltzmann equations 



Here again we only consider leptogenesis from the decays and inverse decays of A'^i. In this approxi- 
mation, the BE for Y/v^ is still given by eq. (1.19) while the BE for 1aLc« lepton asymmetry in the 
flavour a is given b> 



dz 



(1.43) 



Notice that as long as L violation from sphalerons is neglected (see section 1.4) the BE for Yal^ sue 
independent of each other, and hence the solutions for the weak and strong washout regimes are given 
respectively by eq. (1.25)) and eq. (1.27), after replacing ei — >• eia and Ki — )• Kia = Pia^i- 

As an example let us assume that leptogenesis occurs around T ~ 10^° GeV, that is in the two- 
flavour regime. Due to the fast r Yukawa interactions ^i(^i) gets projected onto iriK) ^-nd a coherent 
mixture of e + n eigenstate £e+^ii(^'e+fj}- For illustrative purpose, here we consider a scenario in which 
lepton flavour effects are most prominent. We take both /Ci^,i('ie+^ ^ 1, so that both Ial^ and 
Yal^+u a^re in the strong regime. From eq. (1.27) we can write down the solution: 



Ial(oo) 



>AL,(oo) +yAL, 



vr 



vr^ 
SzfKi 



Y^liO) 



eiY^l{0) + 




-Y^m 



APir , APie+^\ 



uzfKi'^^^"' y Pi 



+ 



pO I 



where in the last line we have used eq. ( 1.40[). If P^ ~ P^^^j^^^ then since AP\t + APie+^j 
second term approximately cancels, and eq. (1.44) reduces to 



(1.44) 
: the 



Ial(oo) 



vr 



3z/i^i 



(1.45) 



We see that the final asymmetry is enhanced by a factor of 2 compared to the unflavoured case. If 



there exists some hierarchy between the flavour projectors, then the second term in eq. (1.44) plays 
an important role and can further enhance the asymmetry. For example we can have P\^^~> P\e+pi. 
while APi^ ^ APie+^. In this case, the second term can dominate over the first term. Finally from 
eq. ( |1.44 ) we also notice that leptogenesis with ei = 0, the so-called purely flavoured leptogenesis 
(PFL^TTF g^j^ indeed proceed [691 ESI ESI EH EH]- In this scenario some symmetry has to be imposed 
to realize the condition ei = 0, as for example an approximate global lepton number U{1)l. In the 
limit of exact U{1)l the active neutrinos will be exactly massless. Instead of the seesaw mechanism, 
the small neutrino masses is explained by U{1)l which is slightly broken by a small parameter 
(the "inverse seesaw") [99j which is technically natural since the Lagrangian exhibits an enhanced 
symmetry when // — )■ [lOOj. In the next section, we will discuss another aspect of flavour effects 
which are in particular crucial for PFL. 



^To study the transition between difTerent flavour regimes (from one to two or from two to tliree flavours), a density 
matrix formalism has to be used [681 1841 194| . 

^''This can also refer to the case where the total CP asymmetry is negligible ei 



0. 
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1.3.4 Lepton flavour equilibration 

Another important effect is lepton flavour equilibration (LFE) [93j. LFE processes violate lepton 
flavour but conserve total lepton number e.g. £aH — )• ipH, and can proceed e.g. via off-shell exchange 
of -/V2,3- In thermal equilibrium, LFE processes can quickly equilibrate the asymmetries generated in 
different flavours. In practice this would be equivalent to a situation where all the flavour projectors 



eq. (1.41) are equal, in which case the flavoured BE eq. (1.43) can be summed up into a single BE: 



^ = eiDi{YN,-Y^l)-PlWiYAL, (1.46) 

where P^^ = 1/2 (1/3) in the two (three) flavours regime. In this case the BE is just like the unflavoured 
case but with a reduced washout which, in the strong washout regime, would result in enhancement of 



a factor of 2 (3) in the two (three) flavours regime (c.f. eq. (1.45)). Clearly, LFE can make PFL with 
ei = impotent |93t l56] . Since LFE N2^3 processes scale as T"^ while the Universe expansion scales 
as T^, in spite of the fact that PFL evades the DI bound, they eventually prevent the possibility of 
lowering too much the leptogenesis scale. A generic study in PFL scenario taking into account LFE 
effects concluded that successful leptogenesis still requires Mi ^10* GeV [97j. A more accurate study 
in the same direction recently carried out in Ref. [98], showed that in fact the leptogenesis scale can 
be lowered down to Mi ~ 10^ GeV. 



1.4 Beyond the Basic Boltzmann Equations 



Within factors of a few, the amount of baryon asymmetry that is generated via leptogenesis in A'^i 
decays is determined essentially by the size of the (flavoured) CP asymmetries and by the rates 
of the (flavoured) washout reactions. However, to obtain more precise results (say, within an 0(1) 
uncertainty) several additional effects must be taken into account, and the formalism must be extended 
well beyond the basic BE discussed in the previous Sections. In the following we review some of the 



most important sources of corrections, namely spectator processes (Section 1.4.1 ), scatterings with top 



quarks and gauge bosons (Section [1.4.2 ), thermal effects (Section 1.4.3), contributions from heavier 



RH neutrinos (Section 1.4.4), and we also discuss the role of quantum corrections evaluated in the 



quantum BE approach (Section 1.4.5). Throughout this review we use integrated BE, i.e. we assume 
kinetic equilibrium for all particle species, and thus we use particles densities instead than particles 
distribution functions. Corrections arising from using non-integrated BE have been studied for example 
in Refs. [101^ I102| 1103^ [TOi] . and are generally subleading. 



1.4.1 Spectator processes 

Reactions that without involving violation oi B — L can still affect the final amount of baryon asymme- 
try are classified as "spectator processes" |105l 1106] . The basic way through which they act is that of 
redistributing the asymmetry generated in the lepton doublets among the other particle species. Since 
the density asymmetries of the lepton doublets are what weights the rates of the washout processes, it 
can be expected that spectator processes would render the washouts less effective and increase the effi- 
ciency of leptogenesis. However, in most cases this is not true: proper inclusion of spectator processes 
implies accounting for all the particle asymmetries, and in particular also for the density asymmetry 
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of the Higgs Yah [.106] . This was omitted in Section 1.2 but in fact has to be added to the density 
asymmetry of the leptons Ya^ in weighting for example washouts from inverse decays. Eq. ( 1.20 ) would 
then become: 

^ = e,Di{YN, - Y^D - 2 {Ym + Y^h) Wi (1.47) 

where the factor of two in front of the washout term counts the leptons and Higgs gauge multiplicity. 
Clearly, in some regimes in which Y^£ and Y/\h are not sufficiently diluted by interacting with other 
particles, this can have the effect of enhancing the washout rates and suppressing the efficiency. 

In the study of spectator processes it is fundamental to specify the range of temperature in which 
leptogenesis occurs. This is because at each specific temperature T, particle reactions must be treated 
in a different way depending if their characteristic time scale r (given by inverse of their thermally 
averaged rates) is [HSl llU7j 

(1) much shorter than the age of the Universe: r <C tu{T); 

(2) much larger than the age of the Universe: r ^> tij{T); 

(3) comparable with the Universe age: T^tu(T). 

Spectator processes belong to the first type of reactions which occur very frequently during one expan- 
sion time. Their effects can be accounted for by imposing on the thermodynamic system the chemical 
equilibrium condition appropriate for each specific reaction, that is J2i f^i = J2f l^F, where /i/ denotes 
the chemical potential of an initial state particle, and fip that of a final state particl^^ The numerical 
values of the parameters that are responsible for these reactions only determine the precise tempera- 
ture T when chemical equilibrium is attained but, apart from this, have no other relevance, and do 
not appear explicitly in the formulation of the problem. Reactions of type (2) cannot have any effect 
on the system, since they basically do not occur. All physical processes are blind to the corresponding 
parameters, that can be set to zero in the effective Lagrangian. In most cases this results in exact 
global symmetries corresponding to conserved charges, and these conservation laws impose constraints 
on the particle chemical potentials. Reactions of type (3) in general violate some symmetry, and thus 
spoil the corresponding conservation conditions, but are not fast enough to enforce chemical equilib- 
rium. These are the only reactions that need to be studied by means of BE, and for which the precise 
value of the parameters that control their rates is of utmost importance. 

A simple case to illustrate how to include spectator processes is the one flavour regime at partic- 
ularly high temperatures (say T ^ lO^^GeV). The Universe expansion is fast implying that except 
for processes induced by the large Yukawa coupling of the top and for gauge interactions, all other 
B — L-conserving reactions fall in class (ii) . Then there are several conserved quantities as for example 
the total number density asymmetries of the RH leptons as well as those of all the quarks except the 
top. Since electroweak sphalerons are also out of equilibrium, B is conserved too (and vanishing, if 
we assume that there is no preexisting asymmetry). B = then translates in the condition: 

2YAQ,+YAt = 0, (1.48) 

where Y/^q.^ is the density asymmetries of one degree of freedom of the top SU{2)l doublet and color 
triplet which, being gauge interactions in equilibrium, is the same for all the six gauge components, 
and Ya* is the density asymmetry of the SU (2)^ singlet top. Hypercharge is always conserved, yielding 

YAQ,+2YAt-YAe + YAH = 0. (1.49) 



The relation between chemical potentials and particle density asymmetries is given in eq. (1.851 
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Finally, in terms of density asymmetries chemical equilibrium for the top- Yukawa related reactions 
/"Qa + I^H = A** translates into 

YAQs + lYAH = YAf (1.50) 

We have three conditions for four density asymmetries, which allows to express the Higgs density 
asymmetry in terms of the density asymmetry of the leptons as Yah = |^A£- Moreover, given that 
only the LH lepton degre es of freedom are populated, we have Yal = 21a£ so that the coefficient 
weighting Wi in eq. (1.47) becomes 2 (Ya^ + Yah) = |Yal and the washout is accordingly stronger. 



With decreasing temperatures, more reactions attain chemical equilibrium, and accounting for 
spectator processes becomes accordingly more complicated. When the temperature drops below T ~ 
10^^ GeV, EW sphalerons are in equilibrium, and baryon number is no more conserved. Then the 



condition eq. (1.48) is no more satisfied and, more importantly, the BE eq. (1.47) is no more valid 
since sphalerons violate also lepton number with in-equilibrium rates. However, sphalerons conserve 
B — L, which is then violated only by slow reactions of type (3), and we should then write down a 
BE for this quantity. Better said, since at T ^ 10^^ GeV all the third generation Yukawa reactions, 
including the ones of the r-lepton, are in equilibrium, the dynamical regime is that of two flavours 
in which the relevant quasi-conserved charges are At- = B/3 — and = B/3 — L^^. The fact 
that only two charges are relevant is because there is always a direction in e-/i space which remains 
decoupled from Ni. The corresponding third charge A' is then exactly conserved, its value can be 



set to zero and the corresponding BE dropped. In this regime, the BE corresponding to eq. (1.47) 
becomes: 

-^^ = eaDi{YN,-Y^l)-2{YAi^+YAH)Wi (q = r, e//) . (1.51) 

To rewrite these equations in a solvable closed form, Yae^ , Yai^^ and Yah must be expressed in 
terms of the two charge densities Ya^ and Ya^^- This can be done by imposing the hypercharge 



conservation condition eq. ( 1.49 ) and the chemical equilibrium conditions that, in addition to eq. ( 1.50 ) 



are appropriate for the temperature regime we are considering. They are [106]: 1. QCD sphaleron 
equilibrium; 2. EW sphaleron equilibrium; 3. 6-quark and r-lepton Yukawa equilibrium. The 'rotation' 
from the particle density asymmetries Yae^ , Yah to the charge densities Ya^ can be expressed in 
terms of the A matrix introduced in [67] Yac^ = A^^^^Ya^ {a, f3 = t, efi) and C-vector Yah = C^Ya^ 
introduced in [BU]. For the present case, with the ordering {e^^T) they are 



It is important to stress that in each temperature regime there are always enough constraints (conser- 
vation laws and chemical equilibrium conditions) to allow to express all the relevant particle density 
asymmetries in terms of the quasi-conserved charges Iaq- This is because each time a conservation 
law has to be dropped (like B conservation above) it gets replaced by a chemical equilibrium condi- 
tion (like EW sphalerons equilibrium), and each time the chemical potential of a new particle species 
becomes relevant, it is precisely because a new reaction involving that particle attains chemical equi- 
librium, enforcing the corresponding condition. As regards the quantitative corrections ascribable to 
spectator processes, several numerical studies have conflrmed that they generally remain below order 
one. Thus, differently from flavour effects, for order of magnitude estimates they can be neglected. 
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Figure 1.2: Diagrams for various 2 o 2 scattering processes: (a) scatterings with the top quarks, 
(6), (c) scatterings with the gauge bosons (^4 = B, Wi with i = 1, 2, 3). 



1.4.2 Scatterings and CP violation in scatterings 

Scattering processes are relevant for the production of the A'"! population, because decay and inverse 
decay rates are suppressed by a time dilation factor oc Mi/T. The Ni = Ni particles can be produced 
by scatterings with the top quark in s-channel /^-exchange qf^ — t- N^a and qi'^ — >• A^£a, by t-channel 
ff-exchange in qi^ Nt'^, ql^ — )• Nf^ and by n-channel //-exchange in t^^a — ^ ^Q, see 



the diagrams labeled (a) in Figure 1.2 Several scattering channels with gauge bosons also contribute 
to the production of A'^i. The corresponding diagrams are labeled {b) and (c) in the same figure. 

When the effects of scatterings in populating the A^i degree of freedom are included, for consistency 
CP violation in scatterings must also be included. In doing so some care has to be put in treating 
properly also all the processes of higher order in the couplings {h^X^, g^X'^, where g' is a gauge coupling) 
with an on-shell intermediate state A^i subtracted out. This can be done by following the procedure 
adopted in Ref. |108j . and we refer to that paper for details. 

In first approximation, the CP asymmetry in scattering processes is the same as in decays and 
inverse decays jTUl llU9j . This result was first found in Refs. |lim [7^ llll| for the case of resonant 
leptogenesis, and was later derived in Ref. [70] for the case of hierarchical Nj. A full calculation of the 
CP asymmetry in scatterings involving the top quark was carried out in Ref. |1U8] . and the validity 
of approximating it with the CP asymmetry in decays was analyzed, finding that the approximation 
is generally good for sufficiently strong RH neutrino hierarchies, e.g. M2/M1 ^ 10. Corrections up to 
several tens of percent can appear around temperatures of order T ~ M2/IO, and can be numerically 
relevant in case of milder hierarchies. 

Regarding the scattering processes with gauge bosons such as N£a AH NH — )• Aia or NA — 
(■aH, their effects in leptogenesis were estimated in Ref. |108j under the assumption that it can also 
be factorized in terms of the decay CP asymmetry. However, with respect to scatterings involving 
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the top quark, there is a significant difference that now box diagrams in wliich the gauge boson is 
attached to a lepton or Higgs in the loop of the vertex- type diagrams are also present, leading to more 
complicated expressions that were explicitly calculated in Ref. |112| . There it was shown that the 
presence of box diagrams implies that for scatterings with gauge bosons the CP asymmetry is different 
from the decay CP asymmetry even for hierarchical RH neutrinos. Still, this difference remains within 
a factor of two |112j so that related effects are in any case not very large. In general, it turns out that 
CP asymmetry in scatterings is more relevant at high temperatures (T > Mi) when the scattering 
rates are larger than the decay rate. Hence, it can be of some relevance to the final value of the 
baryon asymmetry when some of the lepton flavours are weakly washed out, and some memory of the 
asymmetries generated at high temperature is preserved in the final result. 

1.4.3 Thermal corrections 

At the high temperatures at which leptogenesis occurs, the light particles involved in the leptogenesis 
processes are in equilibrium with the hot plasma. Thermal effects give corrections to several ingredients 
in the analysis: (i) coupling constants, (ii) particle propagators (leptons, quarks, gauge bosons and 
the Higgs) and (iii) CP-violating asymmetries, which we briefiy discuss below. A detailed study of 
thermal corrections can be found in Ref. |49| . 

Coupling constants 

Renormalization of gauge and Yukawa couplings in a thermal plasma is studied in Ref. |113j . In 
practice, it is a good approximation to use the zero-temperature renormalization group equations for 
the couplings, with a renormalization scale A ~ 2ttT [l9]. The value A > T is related to the fact that 
the average energy of the colliding particles in the plasma is larger than the temperature. 

The renormalization effects for the neutrino couplings are also well known jll4l I115j . In the non- 
supersymmetric case, to a good approximation these effects can be described by a simple rescaling 
of the low energy neutrino mass matrix m(/x) = r ■ m, where 1.2 ^ r ^ 1.3 for lO^GeV^ /i ^ 
10^^ GeV |49| . and therefore can be accounted for by increasing the values of the neutrino mass 
parameters (for example, m) as measured at low energy by « 20%— 30% (depending on the leptogenesis 
scale). In the supersymmetric case one expects a milder enhancement, but uncertainties related with 
the precise value of the top- Yukawa coupling can be rather large (see Figure 3 in Ref. [IS])- 

Decays and scatterings 

In the thermal plasma, any particle with sizable couplings to the background acquires a thermal mass 
which is proportional to the plasma temperature. Consequently, decay and scattering rates get mod- 
ified. Particle thermal masses have been thoroughly studied in both the SM and the supersymmetric 
SM \in\ Ion 11181 I119| I120j . The singlet neutrinos have no gauge interactions, their Yukawa 
couplings are generally small and, during the relevant era, their bare masses are of the order of the 
temperature or larger. Consequently, to a good approximation, corrections to their masses can be 
neglected. We thus need to account for the thermal masses of the leptons and Higgs doublets and, 
when scatterings are included, also of the third generation quarks and of the gauge bosons (and of 
their superpartners in the supersymmetric case). For a qualitative discussion, it is enough to keep 
in mind that, within the leptogenesis temperature range, we have m//(T) ^ mq^^t{T) ^» m£{T). The 
most important effects are related to four classes of leptogenesis processes: 
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(i) Decays and inverse decays. Since thermal corrections to the Higgs mass are particularly large 
{mH{T) ~ 0.4 r), decays and inverse decays become kinematically forbidden in the temperature range 
in which mniT) — mi{T) < M^-^ < mH{T) + mi{T). For lower temperatures, the usual processes 
A'^i -f-)- IH can occur. For higher temperatures, the Higgs is heavy enough that it can decay: H -f-)- ^A*"!. 
A rough estimate of the kinematically forbidden region yields 2 ^ T/Mi ^ 5. The important point 
is that these corrections are effective only at T > Mi. In the parameter region in > 10~^eV, that is 
favoured by the measurements of the neutrino mass-squared differences, the A^i number density and 
its L-violating reactions attain thermal equilibrium at T ~ Mi and erase quite efficiently any memory 
of the specific conditions at higher temperatures. Consequently, in the strong washout regime, these 
corrections have practically no effect on the final value of the baryon asymmetry. 

(ii) AL = 1 scatterings with top quark. A comparison between the corrected and uncorrected 
rates of the top-quark scattering 7^^^ = ^{q^i ^ iNi) with the Higgs exchanged in the s-channel, 
and of the sum of the t- and u-channel scatterings 7^^^^ = 7(93 -/Vi ^ it) + "fitNi o ^^3) shows 
that the only corrections appearing at low temperatures, and thus more relevant, are for J^f^ (see 
Figure 7.1 in Ref. [109] ). They reduce the scattering rates and suppress the related washouts. This 
peculiar situation arises from the fact that in the zero temperature limit there is a large enhancement 
~ ln(M7Vi/w//) from the quasi-massless Higgs exchanged in the t- and u-channels, which disappears 
when the Higgs thermal mass mH{T) ~ T ~ M^^ is included. 

(iii) AL = 1 scatterings with the gauge bosons. Here the inclusion of thermal masses is required to 
avoid IR divergences that would arise when massless £ (and H) states are exchanged in the t- and u- 
channels. A naive use of some cutoff for the phase space integrals to control the IR divergences can yield 
incorrect estimates of the gauge bosons scattering rates [49j) and would be particularly problematic 
at low temperatures, where gauge bosons scatterings dominate over top-quark scatterings. 

CP asymmetries 

CP asymmetries arise from the interference of tree level and one-loop amplitudes when the couplings 
involved have complex phases and the loop diagrams have an absorptive part. This last requirement 
is satisfied whenever the loop diagram can be cut in such a way that the particles in the cut lines can 
be produced on shell. For the CP asymmetry in decay (at zero temperature) this is guaranteed by 
the fact that the Higgs and the lepton final states coincide with the states circulating in the loops. 
However, in the hot plasma in which A'^i decays occur, the Higgs and the lepton doublets are in thermal 
equilibrium and their interactions with the background introduce in the CP asymmetries a dependence 
on the temperature e — t- e(T) that arises from various effects: 

i ) Absorption and re-emission of the loop particles by the medium require the use of finite temper- 
ature propagators. 

ii) Stimulation of decays into bosons and blocking of decays into fermions in the dense background 
require proper modification of the final states density distributions. 

iii) Thermal motion of the decaying A^'s with respect to the background breaks the Lorentz symmetry 
and affects the evaluation of the CP asymmetries. 

iv) Thermal masses should be included in the finite temperature resummed propagators, and they 
also modify the fermion and boson dispersion relations. Their inclusion yields the most significant 
modifications to the zero temperature results for the CP asymmetries. 
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The first three effects were investigated in Ref. [l8] while the effects of thermal masses was included 
in Ref. |49| . In principle, at finite temperature, there are additional effects related to new cuts that 
involve the heavy A''2,3 neutrino lines. These new cuts appear because the heavy particles in the loops 
may absorb energy from the plasma and go on-shell. However, for hierarchical spectrum, M2,3 ^ Mi, 
the related effects are Boltzmann suppressed by exp(— M2,3/T) that at T ~ Mi is a tiny factor. For 
a non-hierarchical spectrum, the effect of these new cuts can however be sizable. A detail study can 
be found in Ref. [T2T] . 



Propagators and statistical distributions 

Particle propagators at finite temperature are computed in the real time formalism of thermal field the- 
ory [122^1123] . In this formalism, ghost fields dual to each of the physical fields have to be introduced, 
and consequently the thermal propagators have 2x2 matrix structures. For the one-loop computations 
of the absorptive parts of the Feynman diagrams, the relevant propagator components are just those 
of the physical lepton and Higgs fields. The usual zero temperature propagators —iS^{p,me) = 
— mi + ^0+)^-*^ and —iDj^{p, mu) = {p^ — m'jj + iO'^)~^ acquire an additive term that is proportional 
to the particle density distribution rii^H = [ex.p{E£^H /T) ± 

6Se{T) = -2TTne - me) 5{p'^ - mj) , (1.53) 
6DHiT) = +2^ nn 5{p^ - m\) . (1.54) 

For the fermionic thermal propagators, there are other higher order corrections (see Ref. |i49j). Un- 
like the case of bosons, the interactions of the fermions with the thermal bath lead to two different 
types of excitations with different dispersion relations, that are generally referred to as 'particles' and 
'holes' [1^. The contributions of these two fermionic modes were studied in Refs. |124| 11251 1126j where 
it was argued that in the strong washout regime they could give non-negligible effects |126) . The lead- 
ing effects in are proportional to the factor — + nn — 2ni nn that vanishes when the final states 
thermal masses are neglected, because the Bose-Einstein and Fermi-Dirac statistical distributions de- 
pend on the same argument, Ei = Eh = Mi/2. As a consequence, the thermal corrections to the 
fermion and boson propagators {ni and nu) and the product of the two thermal corrections {nenn) 
cancel each other. This was interpreted as a complete compensation between stimulated emission and 
Pauli blocking. As regards the effects in ii), they lead to overall factors that cancel between numerator 
and denominator in the expression for the CP asymmetrj^^ More recently, on the basis of a first 
principle derivation of the CP asymmetry within a quantum BE approach (see Section 1.4.5) it has 
been claimed that the statistical factor induced by thermal loops is instead —ni + nn, which does not 
vanish even in the massless approximation. This would result in a further enhancement in the CP 
asymmetry from the thermal effects |127j . 



Particle motion 

Given that the decaying particle A'^i is moving with respect to the background (with velocity (5) the 
fermionic decay products are preferentially emitted in the direction anti-parallel to the plasma velocity 
(for which the thermal distribution is less occupied), while the bosonic ones are emitted preferentially in 

A similar cancellation holds also in the supersymmetric case. However, because of the presence of the superpartners 
I , H both as final states and in the loops, the cancellation is more subtle and it involves a compensation between the 
two types of corrections i) and ii). We refer to Ref. [48j for details. 
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the forward direction (for which stimulated emission is more effective). Particle motion then induces 
an angular dependence in the decay distribution at order 0{(3). In the total decay rate the 0(/3) 
anisotropy effect is integrated out, and only 0{j3'^) effects remain 0S]. Therefore, while accounting 
for thermal motion does modify the zero temperature results, these corrections are numerically small 
|48[ US] and generally negligible. 

Thermal masses 

When the finite values of the light particle thermal masses are taken into account, the arguments of the 
Bose-Einstein and Fermi-Dirac statistical distributions are different. It is a good approximation [49] 
to use for the particle energies E^^h = Mi/2 =p (m^ — m|)/2Mi. Since now / Eh^ the prefactor 
— n£ + riH — 2n£n}{ that multiplies the thermal corrections does not vanish anymore, and sizable 
corrections become possible. The most relevant effect is that the CP asymmetry vanishes when, as the 
temperature increases, the sum of the light particles thermal masses approaches Mi [l9]. This is not 
surprising, since the particles in the final state coincide with the particles in the loop, and therefore 
when the decay becomes kinematically forbidden, also the particles in the loop cannot go on the mass 
shell. When the temperature is large enough that m//(T) > mg{T) + Mi the Higgs can decay, and then 
there is a new source of lepton number asymmetry associated with H — t- ^A'^i . The CP asymmetry in 
Higgs decays en can be up to one order of magnitude larger than the CP asymmetry in A'^i decays |49] . 
While this could represent a dramatic enhancement of the CP asymmetry, en is non- vanishing only at 
temperatures T ^ Th ~ 5Mi, when the kinematic condition for its decays is satisfied. Therefore, in 
the strong washout regime, no trace of this effect survives. On the other hand, rather large couplings 
are required in order that Higgs decays can occur before the phase space closes: the decay rate can 
attain thermal equilibrium only when m ^ {Tf{/Mi)'^m^ ^ m=K, and therefore, in the weak washout 
regime (m ^ m*), these decays always remain strongly out of equilibrium. This means that only a 
small fraction of the Higgs particles have actually time to decay, and the lepton-asymmetry generated 
in this way is accordingly suppressed. 

In summary, while the corrections to the CP asymmetries can be significant at T ^ Mi (and quite 
large at T ^ Mi for Higgs decays), in the low temperature regime, where the precise value of e plays a 
fundamental role in determining the final value of the baryon asymmetry, there are almost no effects, 
and the zero temperature results still give a reliable approximation. 

1.4.4 Decays of the heavier right-handed neutrinos 

In leptogenesis studies, the effects of A''2,3 are often neglected, which in many cases is not a good 
approximation. This is obvious for example when A^i dynamics is irrelevant for leptogenesis: ei <^ 10^^ 
cannot provide enough CP asymmetry to account for baryogenesis, and rhi <^ m^, implies that A'^i 
washout effects are negligible. It is then clear that any asymmetry generated in decays can 
survive, and becomes crucial for the success of leptogenesis. Another case in which it is intuitively 
clear that effects can be important, is when the RH neutrino spectrum is compact, which means 
that M2,3 have values within a factor of a few from Mi. Then A''i and A''2,3 contributions to leptogenesis 
can be equally important and must be summed up. A model with compact RH neutrino spectrum in 
which A^2,3 dynamics is of crucial importance was recently discussed in Ref. [128j. 

It is less obvious that A''2,3 effects can also be important for a hierarchical RH spectrum and when 
A'^i is strongly coupled. This can happen because decoherence effects related to A'^i-interactions can 
project the asymmetry generated in A''2,3 decays onto a flavour direction that remains protected against 
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A^i washouts \67 \ \129\ \130\ [T3T] . Let us illustrate this with an example. Let us assume that a sizable 
asymmetry is generated in N2 decays, while A^i leptogenesis is inefficient and fails, that is: 

m2 m^,, nil » m*. (1.55) 

Assuming also a strong hierarchy and that leptogenesis occurs thermally guarantees that |131j : 



n7Vi(T~M2) «0, riN^iT Ml) ^ 0. (1.56) 

Thus, the dynamics of and A^i are decoupled: there are neither A'^i-related washout effects during 
leptogenesis, nor A''2-related washout effects during A'"i leptogenesis. The N2 decays into a combination 
of lepton doublets that we denote by £2 ■ 

m = {Y^Y),,'/'J2^M- (1-57) 



The second condition in eq. (1.55) implies that already at T ^ Mi the A^i-Yukawa interactions are 
sufficiently fast to quickly destroy the coherence of £2- Then a statistical mixture of £1 and of the 
state orthogonal to £1 builds up, and it can be described by a suitable diagonal density matrix. Let 
us consider the simple case where both and A^i decay at T ^ 10^^ GeV, so that flavour effects 
are irrelevant. A convenient choice for an orthogonal basis for the lepton doublets is {£i,£o,£q) where, 
without loss of generality, £q satisfies {£o\£2) = 0. Then the asymmetry AYi^ produced in A''2 decays 
decomposes into two components: 

AYe, = AYe„ AY^, = AY,,, (1.58) 

where = |(^o|^2)P and = 1 — c^. The crucial point here is that we expect, in general, 7^ and, 
since (^ol^i) = 0) -^^^o protected against A''i washout. Consequently, a finite part of the asymmetry 
AY^j from decays survives through A^i leptogenesis. A more detailed analysis |131| finds that AY^^ 
is not entirely washed out, resulting in the final lepton asymmetry Yal = (3/2) AY^^ = (3/2)c^ AY^,. 

For 10^ GeV^ Mi ^ 10^^ GeV, flavour issues modify the quantitative details, but the qualitative 
picture, and in particular the survival of a finite part of AY,, , still holds. In contrast, for Mi ^ 10^ GeV, 
the full flavour basis (^ei^^^^r) is resolved and thus there are no directions in flavour space where an 
asymmetry is protected, so that Y^, can be erased entirely. A dedicated study in which the various 
flavour regimes for A'^1,2,3 decays are considered can be found in Ref. [132J. 

In conclusion A'^2,3 leptogenesis cannot be ignored, unless one of the following conditions holds: 

1. The reheat temperature is below M2. 

2. The asymmetries and/or the washout factors vanish, €^2^2 ^ and eN^ils ~ 0- 

3. A''i-related washout is still significant at T ^ 10^ GeV. 



1.4.5 Quantum Boltzmann equations 

So far we have analyzed the leptogenesis dynamics by adopting the classical BE of motion. An 
interesting question which has attracted some attention recently \133\ [85l 1134^ I135| 11361 1127[ 11371 IMl 
I138t I12H [T39J is under which circumstances the classical BE can be safely applied to get reliable results 
and, conversely, when a more rigorous quantum approach is needed. Quantum BE are obtained starting 
from the non-equilibrium quantum field theory based on the Closed Time-Path (CTP) formulation 
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|140j . Both, CP violation from wave- function and vertex corrections are incorporated. Unitarity 
issues are resolved and an accurate account of all quantum-statistical effects on the asymmetry is 
made. Moreover, the formulation in terms of Green functions bears the potential of incorporating 
corrections from Thermal Field Theory within the CTP formalism. 

In the CTP formalism, particle number densities are replaced by Green's functions obeying a set of 
equations which, under some assumptions, can be reduced to a set of kinetic equations describing the 
evolution of the lepton asymmetry and the RH neutrinos. These kinetic equations are non-Mar kovian 
and present memory effects. In other words, differently from the classical approach where every 
scattering in the plasma is independent from the previous one, the particle abundances at a given 
time depend upon the history of the system. The more familiar energy- conserving delta functions are 
replaced by retarded time integrals of time-dependent kernels and cosine functions whose arguments 
are the energy involved in the various processes. Therefore, the non-Markovian kinetic equations 
include the contribution of coherent processes throughout the history of the kernels and the relaxation 
times are expected to be typically longer than the one dictated by the classical approach. 

If the time range of the kernels are shorter than the relaxation time of the particles abundances, 
the solutions to the quantum and the classical BE differ only by terms of the order of the ratio of the 
timescale of the kernel to the relaxation timescale of the distribution. In thermal leptogenesis this is 
typically the case. However, there are situations where this does not happen. For instance, in the case 
of resonant leptogenesis, two RH (s)neutrinos A^i and are almost degenerate in mass and the CP 
asymmetry from the decay of the first RH neutrino A^i is resonantly enhanced if the mass difference 
AM = (M2 — Ml) is of the order of the decay rate of the second RH neutrino r7V2 • The typical 
timescale to build up coherently the CP asymmetry is of the order of 1/AM, which can be larger than 
the timescale ~ ^/^Ni for the change of the abundance of the A^i's. 

Since we need the time evolution of the particle asymmetries with definite initial conditions and not 
simply the transition amplitude of particle reactions, the ordinary equilibrium quantum field theory 
at finite temperature is not the appropriate tool. The most appropriate extension of the field theory 
to deal with non-equilibrium phenomena amounts to generalizing the time contour of integration 
to a closed time-path. More precisely, the time integration contour is deformed to run from to to 
+00 and back to to- The CTP formalism is a powerful Green's function formulation for describing 
non-equilibrium phenomena in field theory. It allows to describe phase-transition phenomena and 
to obtain a self-consistent set of quantum BE. The formalism yields various quantum averages of 
operators evaluated in the in-state without specifying the out-state. On the contrary, the ordinary 
quantum field theory yields quantum averages of the operators evaluated with an in-state at one end 
and an out-state at the other. 

For example, because of the time-contour deformation, the partition function in the in-in formalism 
for a complex scalar field is defined to be 



Z 



= Tr T ^exp 
= Tr 7+ ^exp 
X T- ( exp 



i / (j(x)0(x) + jt(x)(A"^( x)] d X 
i J (j+(x)0+(x) + j|(x)(/>^(x)) d^x 
i j (j_(2;)0_(2;) + j1(x)(/>^x)) d^x 



P 



(1.59) 



where C in the integral denotes that the time integration contour runs from to to plus infinity and 
then back to to again. The symbol p represents the initial density matrix and the fields are in the 
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Heisenberg picture and defined on this closed time-contour (plus and minus subscripts refer to the 
positive and negative directional branches of the time path, respectively). The time-ordering operator 
along the path is the standard one (7+) on the positive branch, and the anti-time-ordering (71) on the 
negative branch. As with the Euclidean-time formulation, scalar (fermionic) fields (p are still periodic 
(anti-periodic) in time, but with (j){t,x) = 4>{t — il3,x), f3 = 1/T. The temperature T appears due to 
boundary condition, and time is now explicitly present in the integration contour. 

We must now identify field variables with arguments on the positive or negative directional branches 
of the time path. This doubling of field variables leads to six different real-time propagators on the 
contour. These six propagators are not independent, but using all of them simplifies the notation. For 
a generic charged scalar field (p they are defined as 

G>(x,y) = -G^+{x,y) = -i{<l>{x)<p\y)), 

G<{x,y) = -G^-{x,y) = -i{cp\y)^{x)), 

G%x,y) = G++{x,y) = 6{x,y)G>{x,y) + 9{y,x)G<{x,y), 

Gl{x,y) = G^-{x,y) = e{y,x)G>{x,y) + 9{x,y)G<{x,y), 

Glf){x, y) = — = G^ — G^, G^(x, y) = G*^ — G^ = G^ — G^, (1.60) 

where the last two Green's functions are the retarded and advanced Green's functions respectively 
and 6{x, y) = 9{tx — ty) is the step function. 

For a generic fermion field ip the six different propagators are analogously defined as 

G^{x,y) = -G;^+ (x,y) = -i(V'(x)^/;(y)), 
G<{x,y) = -G^-{x,y) = +imy)i^{x)), 
G\{x,y) = G++{x,y) = e{x,y)G>{x,y) + e{y,x)G<{x,y), 
G%{x,y) = G;^- ix,y) = e{y,x)G^ix,y) + eix,y)G^{x,y), 

G\i,{x,y) = Glp — G^ = G-^ — G^p, G'^{x,y) = G^p — G^ = G^ — G\p. (1-61) 
From the definitions of the Green's functions, one can see that the hermiticity properties 

(i7°G^(x,y))^ = i70G^(2/,x), (iG^x, y))t = iG^y, x), (1.62) 

are satisfied. For interacting systems, whether in equilibrium or not, one must define and calculate 
self-energy functions. Again, there are six of them: S*, S*, S^, S^, S*" and S". The same relationships 
exist among them as for the Green's functions in (1.60) and (1.61), such as 

E"^ = S* - S< = S> - S*, S" = S* - S> = $]< - S*. (1.63) 

The self-energies are incorporated into the Green's functions through the use of Dyson's equations. A 
useful notation may be introduced which expresses four of the six Green's functions as the elements 
of two-by-two matrices 



G= r.> ^f h S= ( > . , (1.64) 
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where the upper signs refer to the bosonic case and the lower signs to the fermionic case. For systems 
either in equihbrium or in non-equihbrium, Dyson's equation is most easily expressed by using the 
matrix notation 

G{x,y) = G\x,y)+ [ d% / G^x, ^2)^(^2, y), 



(1-65) 

where the superscript "0" on the Green's functions means to use those for noninteracting system. It 
is useful to notice that Dyson's equation can be written in an alternative form, instead of (1.65), with 
G^ on the right in the interaction terms. 



G{x,y) = G^{x,y)+ (fz^ hfz^G{x, z^)T.{z^, Zi)G\zi,y). 



(1.66) 



Eqs. (1.65) and (1.66) are the starting points to derive the quantum BE describing the time evolution 
of the CP-violating particle density asymmetries. 

To proceed, one has to choose a form for the propagators. For a generic fermion ip (and similarly 
for scalars) one may adopt the real-time propagator in the form G'^(k, tx — ty) in terms of the spectral 
function /3^(k, fco) 



G^^i>(k, tx — ty 



-00 27r 

X { [1 - /^.(fcO)] e{tx - ty) - U{ky{ty -tx)}, 



(1.67) 



where f^{k^) represents the fermion distribution function. Again, particles must be substituted by 
quasiparticles, dressed propagators are to be adopted and self-energy corrections to the propagator 
modify the dispersion relations by introducing a finite width T^{k). For a fermion with chiral mass 
m^, one may safely choose 



/9^(k, A:°) = i{/t + m^) 



_(fcO + + ir^)2 - u:l{k) {k^ -ie- ir^)2 - a;2 (fe)_ 



(1.68) 



where a;2 (/c) = k2 + M|(r) and M^(T) is the effective thermal mass of the fermion in the plasma (not 
a chiral mass). Performing the integration over k^ and picking up the poles of the spectral function 
(which is valid for quasi-particles in equilibrium or very close to equilibrium), one gets 



2u;„ 



|(/: + m^) [1 - /^(w^ - iT^)] e~ 



■i{u)^,-ir^){t:c-ty) 



G<{k,ty-tx) = — {(/; + m^)/^(w^ + ir^)e-^('^*-*W--*-) 



+ 7 (/;-m^)7 I - fJuj^-iV^) 



(1.69) 



where k^ = uj^ and denote the distribution function of the fermion particles and antiparticles, 

respectively. The expressions (1.69) are valid for tx — ty > 0. 

The above definitions hold for the lepton doublets (after inserting the chiral LH projector Pl), as 
well as for the Majorana RH neutrinos, for which one has to assume identical particle and antiparticle 
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distribution functions and insert the inverse of the charge conjugation matrix C in the dispersion 
relation. 

To elucidate further the impact of the CTP approach and to see under which conditions one can 
obtain the classical BE from the quantum ones, one may consider the dynamics of the lightest RH 



neutrino A'^i. To find its quantum BE we start from eq. (1.65) for the Green's function of the RH 
neutrino A'^i 



-^Ni z)Gn^ {z, y) + {x, z)G%^ {z, y) 



d^z 



dt. 



(1.70) 



Adopting the corresponding form for the RH neutrino propagator and the center-of-mass coordinates 



X = {t,X) = ]^{x + y), {T,r) = x-y, 



(1.71) 



one ends up with the following equation 



g/iVi(k,t) 
dt 



dt. 



d^ 



P 



1 



1 



/o J (27r)3 2l^^(p) 2ujH{k - p) ujn, (k) 
X [/jv,(k, t)(l - /,(p, t))(l + ///(k - p, t)) 
-Mp,t)fH{k-p,t){l-fN,{k,t))] 
X cos [(wiVj (k) - uji{p) -ujnik- p)) {t - t^)] 



\M{Ni^eH)\' 



~ -2 



dt. 



d^p 



1 



1 



1 



J (27r)3 2oje{p) 2ojH{k - p) u;n, (k) 
x(/^,(k,t)-/;^(p)/^^(k-p)) 
X cos [(wATj (k) - uji{p) -ujnik- p)) {t - t^)] . 



\M{Ni^eH)\' 



(1.72) 



This equations holds under the assumption that the relaxation timescale for the distribution functions 
are longer than the timescale of the non-local kernels so that they can be extracted out of the time 
integral. This allows to think the distributions as functions of the center-of-mass time only. We have 



set to zero the damping rates of the particles in eq. ( 1.69 ) and retained only those cosines giving rise to 



energy delta functions that can be satisfied. Under these assumptions, the distribution function may 
be taken out of the time integral, leading - at large times - to the so-called Markovian description. 



The kinetic equation (1.72) has an obvious interpretation in terms of gain minus loss processes, but the 



retarded time integral and the cosine function replace the familiar energy- conserving delta functions. In 
the second passage, we have also made the usual assumption that all distribution functions are smaller 
than unity and that those of the Higgs and lepton doublets are in equilibrium and much smaller than 
unity, fefn ^ fPf';'' 

any distribution function we may write / 



. Elastic scatterings are typically fast enough to keep kinetic equilibrium. For 



{n/n^'^)f^'^, where n denotes the total number density. 



Therefore, eq. (1.72) can be re- written as 
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dt 



-{rNAt))nN, + {rNAt))nl, 



•eg _ 



COS [(wjvi - - uh) {t - tz)] , 



dt. 



(2^)3 n^^ 



cos [(Wivi - - Wif) - tz)] , 



(1.73) 



where (rjVi(i)) is the time-dependent thermal average of the Lorentz-dilated decay width. Integrating 
over large times, t — >■ oo, thereby replacing the cosines by energy-conserving delta functions 



Jo 



dtz cos [(wiVi -i^t- ujh) {t - tz)] = ttS {uNi -uje- UJh) , 



(1.74) 



we find that the two averaged rates {Tni) and (f at^) coincide and we recover the usual classical BE 
for the RH distribution function 



dt 

d^k f d^p \MiN,^iH)\^ 



' (27r)3 n^^ J [2Tr)-^ 2ui2ujHUJNi 



(1.75) 



Taking the time interval to infinity, namely implementing Fermi's golden rule, results in neglecting 
memory effects, which in turn results only in on-shell processes contributing to the rate equation. The 
main difference between the classical and the quantum BE can be traced to memory effects and to the 
fact that the time evolution of the distribution function is non-Markovian. The memory of the past 
time evolution translates into off-shell processes. 

Similarly, one can show that the equation obeyed by the asymmetry reads 



dt 



-J dtz Tr [S> {X, z)G< {z, X) - G> (X, 2)S< {z, X) 

Gt {X, ^)S> {z, X) - S< {X, z)Gl {z, X)] . (1.76) 



Proceeding as for the RH neutrino equation one finds (including for the moment only the 1-loop wave 
contribution to the CP asymmetry 
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(r 



iVi, 



13=1 



dh 



t2 



-(r.^+rn) 



(i2-ti) 



1 



P)^Afi(k) 



eq. 



d3q 1 - /;^^(q) + /^^(k - q) 
(2^)3 2a;£jq)2wH(k - q)a;7V2(k) 



sin (a;iv-^(t - ti) + (^w^^ + ujh) {h - t2) + ^N2{'t2 - tz) + (wf, + wh) (t^ - 
Tt{MiPl A), 



(1.77) 



where, to avoid double counting, we have not inserted the decay rates in the propagators of the initial 
and final states and, for simplicity, we have assumed that the damping rates of the lepton doublets and 
the Higgs field are constant in time. This should be a good approximation as the damping rate are to 
be computed for momenta of order of the mass of the RH neutrinos. As expected from first principles, 
we find that the CP asymmetry is a function of time and its value at a given instant depends upon 
the previous history of the system. 

Performing the time integrals and retaining only those pieces which eventually give rise to energy- 
conserving delta functions in the Markovian limit, we obtain 

3 

13=1 



\{t) 



dt. 



cos [{loni - '^tc. - ^h) {t - t^)] 



d^k f%\ 

(2^)3 
sin [{uJN., 




p + Tuf + (^Afi 

{^N2 ■ 



WAT J sm 



l-/g(p) + /^^(k-p) 
(27r)3 2w^ Jp)2t^jy(k - p)a;ivi(k) 



d3q 1 _ /;j(q) + /-?(k - q 



Tr(MiPL i>M2 k)- 



(1.78) 



(27r)3 2u;£jq)2wH(k - q)a;^2(k) 

From this expression it is already manifest that the typical timescale for the building up of the coherent 
CP asymmetry depends crucially on the difference in energy of the two RH neutrinos. 

If we now let the upper limit of the time integral to take large values, we neglect the memory 
effects, the CP asymmetry picks contribution only from the on-shell processes. Taking the damping 
rates of the lepton doublets equal for all the flavours and the RH neutrinos nearly at rest with respect 
to the thermal bath, the CP asymmetry reads (now summing over all flavour indices) 

, 2 



it) 



Im {yY'^^ 



12 



Ml 



1 



{YY^\^{YY^\^M2 



N2 



(AM)2 + r^^ 





-AMt' 


[2 AM sin^ 


2 



Tn2 sin [AMt] 



(1.79) 



where AM = {M2 — Mi). The CP asymmetry (1.79) is resonantly enhanced when AM ~ Tj^f^ and at 
the resonance point it is given by 
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1 



(1 - sin [AMt] - cos [AMt]) , 



(1.80) 



The timescale for the building up of the CP asymmetry is ~ 1/AM. The CP asymmetry grows 
starting from a vanishing value and, for t ^ (AM)~^, it averages to the constant standard value. 
This is true if the timescale for the other processes relevant for leptogenesis is larger than ~ 1/AM. 
In other words, one may define an "average" CP asymmetry 



1 



dt'6W(t') 



(1.81) 



where r, 
AL 



the 

the oscillating functions in (|1.80|) are averaged to zero 



p represents the typical timescale of the other processes relevant for leptogenesis, e.g. 
1 scatterings. If Tp > 1/AM 



and the average CP asymmetry is given by the value used in the literature. However, the expression 



(1.79) should be used when Tp <> 1/AM ~ T 



N2- 



The fact that the CP asymmetry is a function of time is particularly relevant in the case in which 
the asymmetry is generated by the decays of two heavy states which are nearly degenerate in mass 
and oscillate into one another with a timescale given by the inverse of the mass difference and has 
been studied in Refs. 



[ I141j . From eq. (1.79) it is manifest that the CP asymmetry itself oscillates 
with the very same timescale and such a dependence may or may not be neglected depending upon the 
rates of the other processes in the plasma. If r^Vi ^ T^^, the time dependence of the CP asymmetry 



may not be neglected. The expression (1.79) can also be used, once it is divided by a factor 2 (because 



in the wave diagram also the charged states of Higgs and lepton doublets may propagate) and the 
limit M2 ^ Ml is taken, for the CP asymmetry contribution from the vertex diagram 



evit) 



Im 



(rr.) 



12 



Ml 



IGvr (W)^^ M2 



2sin^ 



Mot 



Mo 



sin [M2t] 



(1.82) 



In this case, the timescale for this CP asymmetry is ~ M2 and much larger than any other timescale in 
the dynamics. Therefore, one can safely average over many oscillations, getting the expression present 
in the literature. 

What discussed here provides only one example for which a quantum Boltzmann approach is 
needed. In general, the lesson is that quantum BE are relevant when the typical timescale in a quantum 
physical process, such as the timescale for the unflavour-to-flavour transition or the timescale to build 
up coherently the CP asymmetry (of the order of 1/AM) is larger than the timescale for the change 
of the abundances. 



1.5 Supersymmetric Leptogenesis 
1.5.1 What's new? 

Supersymmetric leptogenesis constitutes a theoretically appealing generalization of leptogenesis for 
the following reason: while the SM equipped with the seesaw provides the simplest way to realize 
leptogenesis, such a framework is plagued by an unpleasant fine-tuning problem. For a non degenerate 
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spectrum of heavy Majorana neutrinos, successful leptogenesis requires generically a scale for the 
singlet neutrino masses that is much larger than the EW scale [61] but at the quantum level the 
gap between these two scales becomes unstable. Low-energy supersymmetry (SUSY) can naturally 
stabilize the required hierarchy, and this provides a sound motivation for studying leptogenesis in the 
framework of the supersymmetrized version of the seesaw mechanism. Supersymmetric leptogenesis, 
however, introduces a certain conflict between the gravitino bound on the reheat temperature and the 
thermal production of the heavy singlets neutrinos |142| 11431 11441 . In this section, we will not be 
concerned with the gravitino problem, nor with its possible ways out but focus on the new features 
SUSY brings in for leptogenesis. 

The supersymmetric type-I seesaw model is described by the superpotential of the Minimal Su- 
persymmetric SM (MSSM) with the additional terms: 

w = ]^Mp^n;n'^ + KpN; i^Hu, (i.83) 

where p,q = 1,2,... label the heavy singlet states in order of increasing mass, and a = e, fi,T labels the 



lepton flavour. In eq. (1.83) i, Hu and A^'^, are respectively the chiral superfields for the lepton and the 
up-type Higgs SU{2)i doublets and for the heavy SU{2)l singlet neutrinos defined according to usual 
conventions in terms of their LH Weyl spinor components (for example the N'^ supermultiplet has 
scalar component A^* and fermion component Nf^). Finally the SU{2)l index contraction is defined 
as iaHu = Cpcri^H^ with eu = +1. 

Originally, the issue of MSSM leptogenesis was approached in conjunction with SM leptogenesis [511 
139] as well as in dedicated studies [1461 1147j . However in these first analysis, several features that 
are specific of SUSY in the high temperature regime relevant for leptogenesis, in which soft SUSY 
breaking parameters can be effectively set to zero, had been overlooked. In that case, in spite of the 
large amount of new reactions, the differences between SM and MSSM leptogenesis can be resumed 
by means of simple counting of a few numerical factors [1481 11301 1109] , like for example the number 
of relativistic degrees of freedom in the thermal bath, the number of loop diagrams contributing to 
the CP asymmetries, the multiplicities of the final states in the decays of the heavy neutrinos and 
sneutrinos and one can estimate [109j 

Yab (oo)^^^^ / (strong washout); 



>AB (oo) 



2^/2 (weak washout). 



Recently it was pointed out in Ref. [89] that in fact MSSM leptogenesis is rich of new and non- 
trivial features, and genuinely different from the simpler realization within the SM. The two important 
effects are: 

(a) If the SUSY breaking scale does not exceed by much 1 TeV, above T ~ 5 x 10^ GeV the particle 
and superparticle density asymmetries do not equilibrate |149j . and it is mandatory to account in the 
BE for the differences in the number density asymmetries of the boson and fermion degrees of freedom. 
This can be given in terms of a non-vanishing gaugino chemical potential. 

(b) When soft SUSY breaking parameters are neglected, additional anomalous global symmetries that 
involve both SU{2)l and SU{3)c fermion representations emerge [150]. As a consequence, the EW 
and QCD sphaleron equilibrium conditions are modified with respect to the SM, and this yields a 
different pattern of sphaleron induced lepton-flavour mixing |67 t 1691 170] . In addition, a new anomaly- 
free exactly conserved 7?.-symmetry provides an additional constraint that is not present in the SM and 
a careful counting reveals that four independent quantities, rather than the three of the SM case, are 
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required to give a complete description of the various particle asymmetries in the thermal bath, with 
the additional quantity corresponding to the number density asymmetry of the heavy scalar neutrinos. 

Although the modifications above are interesting from the theoretical point of view, a quantitative 
comparison with the results obtained when the new effects are ignored shows that the corrections 
remain below 0(1) 



89 



13 



Finally, it should be pointed out that in the supersymmetric case, the 



temperature in which the lepton flavour effects (see Section 1.3) come into play is enhanced by a 
factor of (1 + tan^ f3) since the charged Yukawa couplings are given by ha = rria/ {vu cos 

The purpose of the following sections is twofold. We describe the chemical equilibrium conditions 



and conservation laws for MSSM in conjunction to SM. In Section 1.5.2 we list the constraints that 
hold independently of assuming a regime in which particle and sparticle chemical potentials equilibrate 
(superequilibration (SE) regime) or do not equilibrate (non-superequilibration (NSE) regime). In 
Section [l. 5. 3| we list the constraints that hold only in the SE regime, and in Section [l.5.4 the ones that 
hold in NSE regime. The question of NSE is irrelevant in the SM since there are no superparticles. 



Hence the parts relevant for the SM are Section 1.5.2 and Section 1.5.3, with the chemical potential of 
the gaugino set to zero, the chemical potential of the down-type higgsino replaced by the minus of the 



up-type higgsino (see eqs. (1.95) and (1.96)), and all the quantities related to superparticles replaced 
by the ones for particles. 



1.5.2 General constraints 

We first list in items (1), (2) and (3) below the conditions that hold in the whole temperature range 
My/ <C T ^ lO^'^ GeV. Conversely, some of the Yukawa coupling conditions given in items (4) and 
(5) will have to be dropped as the temperature is increased and the corresponding reactions go out of 
equilibrium. First we will relate the number density asymmetry of a particle An = n — n for which a 
chemical potential can be defined to its chemical potential. For both bosons (6) and fermions (/) this 
relation acquires a particularly simple form in the relativistic limit rrihj <C T, and at first order in the 
chemical potential fJ-bj/T <C 1: 



An 



/ 



6 



(1.85) 



For simplicity of notations, in the following we denote the chemical potentials with the same notation 
that labels the corresponding field: (j) = /^c^- 

(1) At scales much higher than Mw, gauge fields have vanishing chemical potential W = B = g = 
[571 ■ This also implies that all the particles belonging to the same SU{2) l or SU (3)c multiplets 
have the same chemical potential. For example (t){Iz = +5) = — ~\) ^ field (j) that is a 
doublet of weak isospin /, and similarly for color. 

(2) Denoting by Wr, Br and (Jr the RH winos, binos and gluinos chemical potentials, and by 
i, Q {i, Q) the chemical potentials of the (s)lepton and (s)quarks LH doublets, the following 
reactions: Q + Qr ^ Q, Q + Wr Q, £ + Wr i, i + Br — £, imply that all gauginos have 
the same chemical potential: —g = Q — Q = —W = i — £ = —B, where W, B and g denote the 
chemical potentials of LH gauginos. It follows that the chemical potentials of the SM particles 



'This modification would be important for certain supersymmetric leptogenesis scenarios which contain new sources 



of CP violation e.g. soft leptogenesis (see Section 1.6.21. 
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are related to those of their respective superpartners as 

Q,l = Q,i + g (1.86) 
Hu,d = Hu,d + 9 (1-87) 
u,d,e = u,d,e — g. (1.88) 

The last relation, in which u,d,e = ur, dn, cr denote the RH SU{2)l singlets, follows e.g. from 
ul = ul + g for the corresponding LH fields, together with = —ur, and from the analogous 
relation for the SU (2) i singlet squarks. 

(3) Before EW symmetry breaking hypercharge is an exactly conserved quantity, and we can assume 
a vanishing total hypercharge for the Universe: 

E (Yaq, + 2^A«, - YAdJ - E (Ym^ + ^AeJ + Y^f,^ - Y^f,^ = 0. (1.89) 

i a 

(4) When the reactions mediated by the lepton Yukawa couplings are faster than the Universe 
expansion ratj^ the following chemical equilibrium conditions are enforced: 

£a-ea + Hd + g = 0, {a = e,fi,T). (1.90) 



If the temperature is not too low lepton flavour equilibration (see Section 1.3.4) induced by 
off-diagonal slepton soft masses will not occur. We assume that this is the case, and thus we 
take the three ia to be independent quantities. 

(5) Reactions mediated by the quark Yukawa couplings enforce the following six chemical equilibrium 
conditions: 

Qi-Ui + Hu + g = 0, {ui = u,c,t), (1-91) 
Qi-di + Hd + g = 0, {di = d, s, b) . (1.92) 

The up-quark Yukawa coupling maintains chemical equilibrium between the LH and RH up-type 
quarks up to T ~ 2 • 10^ GeV. Note that when the Yukawa reactions of at least two families of 
quarks are in equilibrium, the mass basis is fixed for all the quarks and squarks. Intergeneration 
mixing then implies that family-changing charged-current transitions are also in equilibrium: 
bi CL and ti sl imply Q2 = Qs] sl ul and cl — > di imply Qi = Q2. Thus, up to 
temperatures T ^ 10^^ GeV, that are of the order of the equilibration temperature for the charm 
Yukawa coupling, the three quark doublets have the same chemical potential: 

Q = Q3 = Q2 = Qi. (1.93) 

At higher temperatures, when only the third family is in equilibrium, we have instead Q = 
Q3 = Q2 7^ Qi- Above T ~ 10^^ when (for moderate values of tan/3) also 5-quark (as well as 
the r-lepton) SU{2)l singlets decouple from their Yukawa reactions, all intergeneration mixing 
becomes negligible and Q3 Q2 7^ Qi- 



See Section 1.3.1 for the temperature regime when lepton Yukawa interactions are in thermal equilibrium. 
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1.5.3 Superequilibration regime 



At relatively low temperatures, additional conditions from reactions in chemical equilibrium hold. 
Since the constraints below apply only in the SE regime, we number them including this label. 

6sE. Equilibration of the particle-sparticle chemical potentials /x,^ = fi^ is ensured when reactions 

like U —7- U are faster than the Universe expansion rate. These reactions are induced by gaugino 
interactions, but since they require a gaugino chirality flip they turn out to be proportional to 
its mass m^, and can be neglected in the limit rrig — t- 0. 

Furthermore, since the jj, parameter of the HuHd superpotential term is expected to be of the 
order of the soft gaugino masses, it is reasonable to consider in the same temperature range also 
the effect of the higgsino mixing term, which implies that the sum of the up- and down- higgsino 
chemical potentials vanishes. The rates of the corresponding reactions, given approximately by 
m?/T and ~ 1"^ ■, are faster than the Universe expansion rate up to temperatures 



T ^ 5 • 10 



500 GeV 



2/3 



GeV. 



The corresponding chemical equilibrium relations enforce the conditions: 

Hu + iid = 0. 



(1.94) 



(1.95) 
(1.96) 



7sE- Up to temperatures given by eq. (1.94) the MSSM has the same global anomalies than the SM, 



that are the EW SU (2) l-U (1) b+l mixed anomaly and the QCD chiral anomaly. They generate 
the effective operators Oew = ^aiQQQia) [IhXl and Oqcd = Ili{QQu%dli) [T53l fT5T] . 
Above the EW phase transition reactions induced by these operators are in thermal equilibrium, 
and the corresponding conditions read: 



9Q + J2^a = 

a 

6Q-J2(.Ui + di) = 0, 



(1.97) 
(1.98) 



where we have used the same chemical potential for the three quark doublets (eq. (1.93)), which 



is always appropriate in the SE regime below the limit eq. (1.94). 



Eqs. (1.90) and (1.91)-(1.93), together with the SE conditions eqs. (1.95)-(1.96), the two anomaly 
conditions eqs. (1.97)-(1.98) and the hypercharge neutrality condition eq. (1.89), give 11+2+2+1 = 16 



constraints for the 18 chemical potentials. Note however that there is one redundant constraint, that 



we take to be the QCD sphaleron condition, since by summing up eqs. (1.91) and (1.92) and taking 
into account eqs. (1.93), (1.95), and (1.96) we obtain precisely eq. (1.98). Therefore, like in the 
SM, we have three independent chemical potentials. We can define three linear combinations of the 
chemical potentials corresponding to the SU{2)l anomaly-free flavour charges Aq, = B/3 — La that 
being anomaly-free and perturbatively conserved by the low energy MSSM Lagrangian, evolve slowly 
because the corresponding symmetries are violated only by the heavy Majorana neutrino dynamics. 
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Their evolution needs to be computed by means of three independent BE. In terms of the abundances 



eq. (1.18) the density of the charges normalized to the entropy density can be written as: 



1 



(1.99) 



The expression above is completely general and holds in all temperature regimes, including the NSE 



regime (see Section 1.5.4). 



The density asymmetries of the doublet leptons and higgsinos, that weight the washout terms in 
the BE, can now be expressed in terms of the anomaly-free charges by means of the A matrix and C 
vectors introduced respectively in Ref. |i67j and Ref. \69\ that are defined as: 



a 



Ya^. 



(1.100) 



Here and in the following we will give results for the A and C matrices for the fermion states. We 
recall that in the SE regime the density asymmetry of a scalar boson that is in chemical equilibrium 
with its fermionic partner is given simply by y^fe = ^^a/ with the factor of 2 from statistics. 



First generation Yukawa reactions out of equilibrium (SE regime) 

As an example let us now consider the temperatures T ^ 4 • 10^(1 + tan^ /3) GeV, when the d-quark 
Yukawa coupling can be set to zero (in order to remain within the SE regime we assume tan /3 ~ 1). In 
this case the equilibrium dynamics is symmetric under the exchange u -H- d (both chemical potentials 



enter only the QCD sphaleron condition eq. (1.98) with equal weights) and so must be any physical 
solution of the set of constraints. Thus, the first condition in eq. (1.92) can be replaced by the 
condition d = u, and again three independent quantities suffice to determine all the particle density 
asymmetries. The corresponding result is: 



A' 



3 X 2148 



/ 



V 



-906 
75 
75 



120 

-688 
28 




C 



-1 



2148 



(37, 52, 52) 



(1.101) 



Note that since in this regime the chemical potentials for the scalar and fermion degrees of freedom 
of each chiral multiplet equilibrate, the analogous results for Y^£^ + Y^^ can be obtained by simply 
multiplying the A matrix in eq. (1.101) by a factor of 3. This gives the same A matrix obtained 
in the non-supersymmetric case in the same regime (see e.g. eq. (4.12) in Ref. |69j 



The C matrix 

(multiplied by the same factor of 3) differs fro m the non-supersymmetric result by a factor 1/2. This is 
because after substituting = —Hu (see eq. (1.96)) all the chemical potential conditions are formally 
the same than in the SM with Hu identified with the chemical potential of the scalar Higgs, but since 
C expresses the result for number densities, in the SM a factor of 2 from boson statistics appears for 
the SM Higgs. This agrees with the analysis in Ref. [58], and is a general result that holds for SUSY 
within the SE regime. 



1.5.4 Non-superequilibration regime 



At temperatures above the limit given in eq. ( 1.94 ) the Universe expansion is fast enough that reactions 
induced by nig and /i do not occur. Setting to zero in the high temperature effective theory these two 
parameters has the following consequences: 
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Condition eq. (1.95) has to be dropped, and gauginos acquire a non- vanishing chemical potential 
g ^ (corresponding to the difference between the number of LH and RH helicity states). 
The chemical potentials of the members of the same matter supermultiplets are no more equal 



(non-superequilibration) but related as in eqs. (1.86)-(1 



• Condition eq. ( 1.96 ) also has to be dropped, and the chemical potentials of the up- and down-type 
Higgs and higgsinos do not necessarily sum up to zero. 

• The MSSM gains two new global symmetries: rrig — )• yields a global i?-symmetry, while — t- 
corresponds to a global symmetry of the Peccei-Quinn (PQ) type. 

Anomalous and non-anomalous symmetries 

Two linear combinations i?2 and R3 of R and PQ, having respectively only SU{2) l and SU (3)c mixed 
anomalies have been identified in Ref. |150j 

R2 = R-2PQ, Rs = R-3PQ. (1.102) 

The authors of Ref. |15Uj have also constructed the effective multi-fermions operators generated by 
the mixed anomalies: 

Oew = UaiQQQia) HuHdW\ (1.103) 
Oqcd = U.iQQu'd'), f. (1.104) 

Given that three global symmetries B, L and R2 have mixed SU{2)i anomalies (but are free of 
SU{'i)c anomalies) we can construct two anomaly-free combinations, the first one being B — L which 
is only violated perturbatively by N'^ £Hu and the second anomaly- free combination which is also an 
exact symmetry of the MSSM+seesaw in the NSE regime [89| 

TZ=^B-L + R2, (1.105) 

and is exactly conserved. In the SU {3)c sector, besides the chiral anomaly we now have also R3 mixed 
anomalies. Thus also in this case anomaly-free combinations can be constructed, and in particular 
we can define one combination for each quark superfield. Assigning to the LH supermultiplets chiral 
charge x = ~1 these combinations have the form |89j : 

XqL + l^qL ^3, (1.106) 

where, for example, k^c = k^^c =1/3 and kq^ = 2/3. Note that since R3 is perturbatively conserved 
by the complete MSSM-l-seesaw Lagrangian, when the Yukawa coupling of one quark is set to zero the 



corresponding charge eq. (1.106) will be exactly conserved. 



Constraints in the non-superequilibration regime 

In the NSE regime, the conditions listed in items 6se and 7se of the previous section have to be 
dropped, but new conditions arise. 
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^NSE- The conservation law for the TZ charge yields the following global neutrality condition: 

f b 
6 



[ ^ {2Qi - 5ui + Adi) +2j2iia + ea) + 5Hd -Hu + Slg]- An^^ = 0. (1. 



107) 



The last terms in both lines of eq. (1.107) correspond to the contribution to 7^- neutrality from 

-7^^rc = -1. Note 



the lightest sneutrino asymmetry An 



n 



rij^t with charge TZ 



that since in general A'^i is not in chemical equilibrium, no chemical potential can be associated 
to it, and hence this constraint needs to be formulated in terms of its number density asymmetry 
that has to be evaluated by solving a BE for Ya^ = ^^Vi ~ ^N* (^^^ Section 1.5.5). 



Tnse- The operators in eqs. ( 1.103 )-( 1.104) induce transitions that in the NSE regime are in chemical 
equilibrium. This enforces the generalized EW and QCD sphaleron equilibrium conditions |150) : 

3EQ* + E^" + ^« + ^rf + 4^ = 0, (1.108) 



(1.109) 



that replace eqs. (1.97) and (1.98). 



8nse- The chiral-i?3 charges in eq. (1.106) are anomaly-free, but clearly they are not conserved by the 

lark supermultiplet decouples from its Yukawa 
The conservation laws corresponding to these 



quark Yukawa interactions. However, when a quark supermultiplet decouples from its Yukawa 
interactions an exact conservation law arisen 



15 



symmetries read: 



r2 1 

— [3(?R + 6{qR -g)] + - Rs tot 

D 3 

^ 2 [3Ql + 6(Ql+5)] - ^^3 tot 



(1.110) 

6-— 3 ° 

and hold for = Ui, di and Ql = Qi in the regimes when the appropriate Yukawa reactions 
are negligible. Note the factor of 2 for the Ql chiral charge in front of the first square bracket in 
eq. (1.111) that is due to SU{2)l gauge multiplicity. In terms of chemical potentials and An^^, 
the total R3 charge in eqs. ( 1.110| ) and (1.111) reads: 

Rstot = ? (825-3E(2Q* + llM*-4di)+E(16^a + 13e„) + 16#rf-14F„. 



(1.112) 



where R 



-1. As regards the leptons, since they do not couple to the QCD anomaly, by 



setting /le — )• a symmetry under chiral supermultiplet rotations is directly gained for the RH 
leptons implying An^ + Ahe = and giving the condition: 



e q 

3^ 



0. 



(1.113) 



No analogous condition arises for the lepton doublets relevant for leptogenesis, since by assump- 
tion they remain coupled via Yukawa couplings to the heavy A^'s. 

^^Note that hu,d — > implies it and d decoupling, but Qi decoupling is ensured only if also hc,s — > 0. 



36 



Neutrino Physics 



Leptogenesis in the Universe 



In the NSE regime there are different flavour mixing matrices for the scalar and fermion components 
of the leptons and Higgs supermultiplets. To express more concisely all the results, it is convenient to 
introduce a new C vector to describe the gaugino number density asymmetry per degree of freedom 
in terms of the relevant charges: 



with 



(1.114) 



First generation Yukawa reactions out of equilibrium (NSE regime) 



As an example, in the temperature range between 10*^ and 10^^ GeV, and for moderate values of tan /3 



all the first generation Yukawa couplings can be set to zero. Using for n, d conditions eq. (1.110) and 
for e condition eq. (1.113) as are implied by /i^^d, /ig — we obtain: 

/ -198117 



A' 



C9 



9 X 162332 



V 



26634 
26634 



33987 
-147571 
14761 



33987 
14761 
-147571 



-8253 
-8055 
-8055 



-11 



C 



c 



H, 



162332 

-1 
162332 
1 



(163, 165, 165, -255), 
(3918, 4713, 4713, 95) , 

(5413, 9712, 9712, -252) 



(1.115) 



For completeness, in eq. (1.115) we have 



3 X 162332 

where the rows correspond to (Ya^, Ya^, Ya^, Ya 

also given the results for C^'^ even if only the up-type Higgs density asymmetry is relevant for the 
leptogenesis processes. Note that neglecting the contribution of An^^ to the global charges T^tot in 



eq. (1.107) and -Rstot in eq. (1.112) corresponds precisely to setting to zero the fourth column in all the 



previous matrices. Then, analogously with the SE and SM cases, within this '3-column approximation' 
all particle density asymmetries can be expressed just in terms of the three Ya„ charge densities. 

1.5.5 Supersymmetric Boltzmann equations 

In order to illustrate how the new effects described above modify the structure of the BE, here we 



write down a simpler expressions in which only decays and inverse decays are includec 



sHz- 



dz 



sHz- 



sHz 



sHz 



dY, 



dz 

dYA- 
dz 

dYA^ 
dz 



Yn^ 
'Y 



1 7iVi: 



^req 1 Ni n 'A^i ^req ■ • • ; 



(1.116) 

(1.117) 
(1.118) 



^ - 1 j IN, + 




2 hiv, 



The complete set of BE including decays, inverse decays and scatterings with top-quark is given in the Appendix of 



Ref. 
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+Pla 7^1 + 



1 



a 



i 

aa 



(1.119) 



where 7jy^ is the corresponding thermally averaged decay rate for RH sneutrino A'^i. In eq. 



we have introduced the overall sneutrino abundance 



+ while Y^. 



1.117) 



m 



eq. (1.118) is the sneutrino density asymmetry that was already introduced in Section 1.5.4 



washout terms we have normalized the charge densities 
of a fermion with one degree of freedom Yf'^. 



In the 

to the equilibrium density 



In eqs. (1.116) and (1.119) we have also neglected for 



simplicity all finite temperature effects. Taking these effects into account would imply for example 
that the CP asymmetry for N decays into fermions is different from the one for decays into scalars 



while we describe both CP asymmetries with e^. A few remarks regarding eq. (1.118) are in order. 
In the SE regime g = and thus it would seem that the sneutrino density asymmetry Ya^ vanishes. 
However, this only happens for decays and inverse decays, and it is no more true when additional terms 
related to scattering processes, that are represented in the equation by the dots, are also included (see 
Ref. |147j and the Appendix of Ref. [89]). Therefore, also in the SE regime Yj^^ and Y^, in general 
differ. However, in this case recasting their equations in terms of two equations for Y^^ and Y^^ is 
just a convenient parametrization. On the contrary, in the NSE regime this is mandatory, because 
the sneutrinos carry a globally conserved 7^-charge and 1a - is required to formulate properly the 
corresponding conservation law. As we have seen, this eventually results in Ya - contributing to the 
expressions of the lepton flavour density asymmetries in terms of slowly varying quantities. 

In Ref. [89j a complete numerical analysis was carried out and it was shown that numerical correc- 
tions with respect to the case when NSE effects are neglected remain at the 0(1) level. This is because 
only spectator processes get affected, while the overall amount of CP violation driving leptogenesis 
remains the same than in previous treatments. 



1.6 Beyond Type-I Seesaw and Beyond the Seesaw 

There exist many variants of leptogenesis models beyond the standard type-I seesaw. In this section, 
we try to classify them into appropriate groups. Unavoidably there would be some overlap i.e. a 
hybrid model of leptogenesis which can belong to more than one group, e.g. soft leptogenesis (Section 
from resonantly enhanced CP asymmetry could rightly fall under resonant leptogenesis (Section 
However we try our best to categorize them according to the main features of the model and. 



1.6.2 



1.6.1 



when appropriate, they will be quoted in more than one place. Clearly, the number of beyond type-I 
seesaw leptogenesis models is quite large. We have not attempted in any way to be exhaustive, and 
we apologize in advance for the unavoidable several omissions. 

1.6.1 Resonant leptogenesis 

A resonant enhancement of the CP asymmetry in A^i decay occurs when the mass difference between 
A'^i and is of the order of the decay widths. Such a scenario has been termed 'resonant leptogenesis', 
and has benefited from many studies in different formalism^ (SUl E21 ttSHl ESZl [ESI [ESI 
[ml fTM ITHTl fT62l fT63l [ESI [IH [EH MB [Ml [Ml [leg (see Ref. [I69] for a recent review). The 



See Ref. |154j for a comparison of the different calculations. 
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resonant effect is related to the self energy contribution to the CP asymmetry. Consider, for simplicity, 
the case where only is quasi-degenerate with A^^i. Then, the self-energy contribution involving the 
intermediate to the total CP asymmetry (we neglect important flavour effects [7H]) is given by 



..(self - energy) - M^M^ - M?) Im[(AtA)f,j 

6i(seh energy) - {M^ - M^y + Mlr%^ i\n)u{X^ X)22- ^ ' 



The resonance condition reads 



In this case 



\M2-Mi\ = ^. (1.121) 



I ei (resonance) I ~ - tttTV^TTTTV^ • (1.122) 



1 |Xm[(AtA)f2 
2(AtA)ii(AtA)22 

Thus, in the resonant case, the asymmetry is suppressed by neither the smallness of the light neutrino 
masses, nor the smallness of their mass splitting, nor small ratios between the singlet neutrino masses. 
Actually, the CP asymmetry could be of order one (more accurately, |ei| < 1/2). 

With resonant leptogenesis, the BE are different. The densities of A'^i and are followed, since 
both contribute to the asymmetry, and the relevant timescales are different. For instance, the typical 
time-scale to build up coherently the CP asymmetry is unusually long, of order 1/ AM. In particular, it 
can be larger than the time-scale for the change of the abundance of the sterile neutrinos. This situation 
implies that for resonant leptogenesis quantum effects in the BE can be significant [851 l86l \170\ I141j 



(see Section 1.4.5). 



The fact that the asymmetry could be large, independently of the singlet neutrino masses, opens 
up the possibility of low scale resonant leptogenesis. Models along these lines have been constructed in 
Refs. |156t [75 t[16Hll71j . It is a theoretical challenge to construct models where a mass splitting as small 
as the decay width is naturally achieved. For attempts that utilize approximate flavour symmetries 
see, for example, Refs. |157l 1158^ 1160^ 1164^ I172| I166j . while studies of this issue in the framework of 
minimal flavour violation can be found in Refs. |162l I163| 18 7j . The possibility of observing resonant 
CP violation due to heavy RH neutrinos at the LHC was studied in Refs. \l7'i\ 1174] . 

1.6.2 Soft leptogenesis 



The modifications to standard type-I leptogenesis due to SUSY have been discussed in Section 1.5 



The important parameters there are the Yukawa couplings and the singlet neutrino parameters, which 



appear in the superpotential eq. (1.83). SUSY must, however, be broken. In the framework of the 
MSSM extended to include singlet neutrinos (MSSM+N), there are, in addition to the soft SUSY 
breaking terms of the MSSM, terms that involve the singlet sneutrinos Ni, in particular bilinear (B) 
and trilinear (A) scalar couplings. These terms provide additional sources of lepton number violation 
and of CP violation. Scenarios where these terms play a dominant role in leptogenesis have been 
termed 'soft leptogenesis' [I75l ETSl ETTl UZHl [Hal [THOl [Ml EM [Ml [Ml [Ml [Ml [Ml (Ml [Ml 
nm [Ton nm nm nm mm nm nm \W7\ (see also Ref . [SH for a recent review) . 

Soft leptogenesis can take place even with a single RH neutrino because the presence of the B 
term implies that N and states mix to form two mass eigenstates with mass splitting proportional 
to B itself. Furthermore when B ~ F-, the CP asymmetry is resonantly enhanced realizing the 



resonant leptogenesis scenario (see Section 1.6.1). In the following we will consider a single generation 
MSSM+N. The relevant soft SUSY terms involving N, the SU{2)l gauginos A^'^ the U{1)y gauginos 
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Ai and the three sleptons ia in the basis in which charged lepton Yukawa couplings are diagonal are 
given by 

1 

2' 



£soft = M^N*N + ( AY^eabN^Ht + ^BMNN + h.c. 



+1 (m2A^'°PLA^'° + miAiPiAi + h.c.) , (1.123) 

where for simphcity, proportionahty of the bihnear and trihnear soft breaking terms to the corre- 
sponding SUSY invariant couplings has been assumed: Bm = BM and = AYa- The Lagrangian 



derived from eqs. (1.83) and (1.123) is characterized by only three independent physical phases: 
(pA = arg(^i?*), 4>g^ = I arg (57712) a-'^d (pgy = ^arg(i?m*) which can be assigned to A, and to 
the gaugino coupling operators g2, gy respectively. As mentioned earlier, a crucial role in soft lepto- 
genesis is played by the iV — iV^ mixing to form the mass eigenstates 

iV+ = (e**/2^ + e-'*/2^*) , (1.124) 

iV_ = -^(e^*/2jV-e-**/2^*), (1.125) 
v2 

where $ = arg {BM) and the corresponding mass eigenvalues are = + it |i?M|. Without 
loss of generality, we can set <I> = 0, which is equivalent to assigning the phases only to A and Y^. 

It has been pointed out that the CP asymmetries for the decays of A'^-i- into scalars and fermions 
have opposite sign and cancel each other at the leading order |176| I177|, I192j , resulting in a strongly 
suppressed total CP asymmetry ~ 0{m^^j^/M^) where rrisoft is the scale of soft SUSY breaking terms. 
There are two possibilities that can rescue leptogenesis: Firstly, thermal effects which break SUSY 



can spoil this cancellation |176l 11771 1192| . Secondly, non-superequilibration effects (see Section 1.5.4) 
which imply that lepton and slepton asymmetries differ, can also spoil this cancellation [107J. The CP 
asymmetries for the decays of N± into scalars and fermions are respectively given by 

e%{T) = e-„A^(r), (1.126) 
el{T) = -6-„A/(r), (1.127) 

where is the temperature independent term of ~ 0{msoft/M) which contains contributions from 
the self-energy correction, vertex correction and the interference between the two. In the limit T — )• 0, 
we have A*(T), A-^ (T) — )• 1/2, and thus the inclusion of thermal effects and/or non-superequilibration 
is mandatory to avoid the cancellation between the asymmetries into scalars and fermions. 

We can make a rough estimate of the scale relevant for soft leptogenesis by requiring |e| ~ 
rUsoft/M Z 10-6 which gives M <, 10^ GeV for m-so/t ~ 1 TeV. Hence soft leptogenesis always 
happens in the temperature regime where lepton flavour effects are relevant [190j . In general, the 
CP asymmetry from self-energy contribution requires B <^ w-so/t to be resonantly enhanced. How- 
ever it was shown in Ref. [195] that flavour effects can greatly enhance the efficiency and eventually 
B ~ rrigoft is allowed. The nice feature of soft leptogenesis is that the tension with the gravitino 
problem gets generically relaxed and, in the lower temperature window, is completely avoided. 

1.6.3 Dirac leptogenesis 

The extension of the SM with singlet neutrinos allows for two different ways for generating tiny 
neutrino masses. The first one is the seesaw mechanism which has at least three attractive features: 
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• No extra symmetries (and, in particular, no global symmetries) have to be imposed. 

• The extreme lightness of neutrino masses is linked to the existence of a high scale of new physics, 
which is well motivated for various other reasons {e.g. gauge unification). 

• Lepton number is violated, which opens the way to leptogenesis. 

The second way is to impose lepton number and give to the neutrinos Dirac masses. A priori, one 
might think that all three attractive features of the seesaw mechanism are lost. Indeed, one must 
usually impose additional symmetries. But one can still construct natural models where the tiny 
Yukawa couplings that are necessary for small Dirac masses are related to a small breaking of a 
symmetry. What is perhaps most surprising is the fact that leptogenesis could proceed successfully 
even if neutrinos are Dirac particles, and lepton number is not (perturbatively) broken \W7\ 1198] . 
Such scenarios have been termed 'Dirac leptogenesis' [Ml EMI [2001 [Ml W2\ [2031 1201 [2051 [Ml [M] • 

An implementation of the idea is the following. A CP-violating decay of a heavy particle can result 
in a non-zero lepton number for LH particles, and an equal and opposite non-zero lepton number for 
RH particles, so that the total lepton number is zero. For the charged fermions of the SM, the Yukawa 
interactions are fast enough that they quickly equilibrate the LH and the RH particles, and the lepton 
number stored in each chirality goes to zero. This is not true, however, for Dirac neutrinos. The size of 
their Yukawa couplings is A ^ 10^^^, which means that equilibrium between the lepton numbers stored 
in LH and RH neutrinos will not be reached until the temperature falls well below the electroweak 
breaking scale. To see this, note that the rate of the Yukawa interactions, given by F^^ ~ A^T, becomes 
significant when it equals the expansion rate of the Universe, H T^/rripi. Thus, the temperature 
of equilibration between LH and RH neutrinos is T ~ A^rripi ~ (A/10^^^)^ MeV, that is well below 
the temperature when sphalerons, after having converted part of the LH lepton asymmetry into a net 
baryon asymmetry, are switched off. 

A specific example of a supersymmetric model where Dirac neutrinos arise naturally is presented in 
Ref. [199j . The Majorana masses of the iV-superfields are forbidden by U{1)l. The neutrino Yukawa 
couplings are forbidden by a U{l)]\f symmetry where, among all the MSSM+N fields, only the N 
superfields are charged. The symmetry is spontaneously broken by the vacuum expectation value of a 
scalar field x that can naturally be at the weak scale, (x) ~ Vu- This breaking is communicated to the 
MSSM+N via extra, vector-like lepton doublet fields, (j) + (p, that have masses much larger than 
Vu. Consequently, the neutrino Yukawa couplings are suppressed by the small ratio {x)/M^. The CP 
violation arises in the decays of the vector-like leptons, whereby T{(p — )• NH^) ^ T{(j) — )• N'^Hy) and 
T{(f) — >• Lx) 7^ F(0 — )• L^x'^)- The resulting asymmetries in N and in L are equal in magnitude and 
opposite in sign. Finally note that the main phenomenological implication of Dirac leptogenesis is the 
absence of any signal in neutrinoless double beta decays. 

1.6.4 Triplet scalar (type-II) leptogenesis 

One can generate seesaw masses for the light neutrinos by tree-level exchange of SU{2) j^-triplet scalars 
T im [2081 EQU [2101 [211]. The relevant new terms in the Lagrangian are 

Ct = -M^\T\^ + ^ {[XlU ^ai(3T + MtA^ #T* + h.c.) . (1.128) 

Here, Mt is a real mass parameter, A^, is a symmetric 3x3 matrix of dimensionless, complex Yukawa 
couplings, and A^ is a dimensionless complex coupling. Since this mechanism necessarily involves 
lepton number violation and allows for new CP- violating phases, it is interesting to examine it in the 
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light of leptogenesis [2l2l EUl EH EEl EH M EIZl Eig EH UMl E2Ql EZD EZa ESI EMI EH ESSl 

[2261 12271 [2281 [2291 [2301 [2311 [2321 [233] . One obvious problem in this scenario is that, unlike singlet 
fermions, the triplet scalars have gauge interactions that keep them close to thermal equilibrium at 
temperatures T ^ 10^^ GeV. It turns out, however, that successful leptogenesis is possible even at a 
much lower temperature. This subsection is based in large part on Ref. [220j where further details 
and, in particular, an explicit presentation of the relevant BE can be found. 

The CP asymmetry that is induced by the triplet scalar decays is defined as follows: 

i T + i f 

where the overall factor of 2 comes because the triplet scalar decay produces two (anti)leptons. 



To calculate ey, one should use the Lagrangian in eq. (1.128). While a single triplet is enough 
to produce three light massive neutrinos, there is a problem in leptogenesis if indeed this is the only 
source of neutrinos masses: The asymmetry is generated only at higher loops and in unacceptably 
small. It is still possible to produce the required lepton asymmetry from a single triplet scalar decays 
if there are additional sources for neutrino masses, such as type I, type III, or type II contributions 
from additional triplet scalars. Define mji (mi) as the part of the light neutrino mass matrix that 
comes (does not come) from the contributions of the triplet scalar responsible for ct- 

m = mu + mi. (1.130) 

Then, assuming that the particles exchanged to produce mi are all heavier than T, we obtain the the 
CP asymmetry 

^ ^ ^Ag^ Im[Tr(mi\mi)] 
^ 47r t;2 ^ " Tr(mj\mii) 

where Bl (Bh) is the tree-level branching ratio to leptons (Higgs doublets). If these are the only 
decay modes, i.e. Bl + Bh = 1, then Bl/Bh = / [XhX^h) , and there is an upper bound on 

the asymmetry: 

l^^l < I^^^BlBuY.^- (1-132) 

Note that, unlike the singlet fermion case, le^l increases with larger rriy.. 

As concerns the efficiency factor, it can be close to maximal, r/ ~ 1, in spite of the fact that the 
gauge interactions tend to maintain the triplet abundance very close to thermal equilibrium. There 
are two necessary conditions that have to be fulfilled by the decay rates T ^ II and T — )■ (pep in order 
that this will happen [220]: 

1. One of the two decay rates is faster than the TT annihilation rate. 

2. The other decay mode is slower than the expansion rate of the Universe. 

The first condition guarantees that gauge scatterings are ineffective: the triplets decay before anni- 
hilating. The second condition guarantees that the fast decays do not washout strongly the lepton 
asymmetry: lepton number is violated only by the simultaneous presence of T ^ ii and T — )• (p<p. 



Combining a calculation of rj with the upper bound on the CP asymmetry (1.132), successful 
leptogenesis implies a lower bound on the triplet mass Mt varying between 10^ GeV and 10^^ GeV, 
depending on the relative weight of nin and nii in the light neutrino mass. 
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Interestingly, in the supersymmetric framework, "soft leptogenesis" (see Section 1.6.2) can be 
successful even with the minimal set of extra fields - a single T + T - that generates both neutrino 
masses and the lepton asymmetry [18^ I185| . 

1.6.5 Triplet fermion (type-Ill) leptogenesis 

One can generate neutrino masses by the tree level exchange of 5'C/(2)i-triplet fermions T-^ \234:\ 12351 
1236] (z denotes a heavy mass eigenstate while a is an SU{2)i index) with the Lagrangian 

Ct-^ = [XrUr^J'^rn - Im,T,^T,- + h.c. (1.133) 

Here r"^ are the Pauli matrices, Mj are real mass parameters and At is a 3 x 3 matrix of complex 
Yukawa couplings. 

This mechanism necessarily involves lepton number violation, and allows for new CP-violating 
phases so we should examine it as a possible source of leptogenesis [2371 12381 12281 12391 12401 12411 12421 
1243] . This subsection is based in large part on Ref. \238\ where further details and the relevant BE 
can be found. 

As concerns neutrino masses, all the qualitative features are very similar to the singlet fermion case. 
As concerns leptogenesis there are, however, qualitative and quantitative differences. With regard to 
the CP asymmetry from the lightest triplet fermion decay, the relative sign between the vertex and 
self-energy loop contributions is opposite to that of the singlet fermion case. Consequently, in the 
limit of strong hierarchy in the heavy fermion masses, the asymmetry in triplet decays is three times 
smaller than in the decays of the singlets. On the other hand, since the triplet has three components, 
the ratio between the final baryon asymmetry and er] is three times bigger. The decay rate of the 
heavy fermion is the same in both cases. This, however, means that the thermally averaged decay 
rate is three times bigger for the triplet, as is the on-shell part of the AL = 2 scattering rate. 

A significant qualitative difference arises from the fact that the triplet has gauge interactions. The 
effect on the washout factor rj is particularly significant for m <C 10~^ eV, the so-called "weak washout 
regime" (note that this name is inappropriate for triplet fermions). The gauge interactions still drive 
the triplet abundance close to thermal equilibrium. A relic fraction of the triplet fermions survives. 
The decays of these relic triplets produce a baryon asymmetry, with 

rj ^ Mi/10^^ GeV (for m < 10"^ eV). (1.134) 

The strong dependence on Mi results from the fact that the expansion rate of the Universe is slower 
at lower temperatures. On the other hand, for fh ^ 10~^ eV, the Yukawa interactions keep the heavy 
fermion abundance close to thermal equilibrium, so the difference in r] between the singlet and triplet 
case is only 0(1). Ignoring flavour effects, and assuming strong hierarchy between the heavy fermions, 
Ref. |216] obtained the lower bound 

Ml ^ 1.5 X 10^° GeV. (1.135) 

When the triplet fermion scenario is incorporated in a supersymmetric framework, and the soft break- 
ing terms do not play a significant role, the modifications to the above analysis is by factors of 0(1). 
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1.7 Conclusions 

During the last few decades, a large set of experiments involving solar, atmospheric, reactor and accel- 
erator neutrinos have converged to establish that the neutrinos are massive. The seesaw mechanism 
extends the Standard Model in a way that allows neutrino masses, and it provides a nice explanation 
of the suppression of the neutrino masses with respect to the electroweak breaking scale. Furthermore, 
without any addition or modification, it can also account for the observed baryon asymmetry of the 
Universe. The possibility of giving an explanation of two apparently unrelated experimental facts - 
neutrino masses and the baryon asymmetry - within a single framework that is a natural extension 
of the Standard Model, together with the remarkable 'coincidence' that the same neutrino mass scale 
suggested by neutrino oscillation data is also optimal for leptogenesis, make the idea that baryogenesis 
occurs through leptogenesis a very attractive one. 

Leptogenesis can be quantitatively successful without any fine-tuning of the seesaw parameters. 
Yet, in the non-supersymmetric seesaw framework, a fine-tuning problem arises due to the large cor- 
rections to the mass-squared parameter of the Higgs potential that are proportional to the heavy 
Majorana neutrino masses. Supersymmetry can cure this problem, avoiding the necessity of fine tun- 
ing; however, it brings in the gravitino problem |144| that requires a low reheat temperature after 
inflation, in conflict with generic leptogenesis models. Thus, constructing a fully satisfactory theo- 
retical framework that implements leptogenesis within the seesaw framework is not a straightforward 
task. 

From the experimental side, the obvious question to ask is if it is possible to test whether the 
baryon asymmetry has been really produced through leptogenesis. Unfortunately it seems impossible 
that any direct test can be performed. To establish leptogenesis experimentally, we need to produce 
the heavy Majorana neutrinos and measure the CP asymmetry in their decays. However, in the most 
natural seesaw scenarios, these states are simply too heavy to be produced, while if they are light, 
then their Yukawa couplings must be very tiny, again preventing any chance of direct measurements. 

Lacking the possibility of a direct proof, experiments can still provide circumstantial evidence 
in support of leptogenesis by establishing that (some of) the Sakharov conditions for leptogenesis 
are realized in nature. Planned neutrinoless double beta decay (Oz^/3/3) experiments (GERDA |244) . 
MAJORANA [245J, CUORE [MS]) aim at a sensitivity to the effective Ouf3f3 neutrino mass in the 
few X 10 meV range. If they succeed in establishing the Majorana nature of the light neutrinos, 
this will strengthen our confidence that the seesaw mechanism is at the origin of the neutrino masses 
and, most importantly, will establish that the first Sakharov condition for the dynamical generation 
of a lepton asymmetry (L violation) is realized. Proposed SuperBeam facilities |247l I248j and second 
generation off-axis SuperBeam experiments (T2HK [249j . NOz^A |250j ) can discover CP violation in 
the leptonic sector. These experiments can only probe the Dirac phase of the neutrino mixing matrix. 
They cannot probe the Majorana low energy or the high energy phases, but the important point is 
that they can establish that the second Sakharov condition for the dynamical generation of a lepton 
asymmetry is satisfied. As regards the third condition, that is that the heavy neutrino decays occurred 
out of thermal equilibrium, it might seem the most difficult one to test experimentally. In reality the 
opposite is true, and in fact we already know that an absolute neutrino mass scale of the order of the 
solar or atmospheric mass differences is perfectly compatible with sufficiently out of equilibrium heavy 
neutrinos decays. 

Given that we do not know how to prove that leptogenesis is the correct theory, we might ask if there 
is any chance to falsify it. Indeed, future neutrino experiments could weaken the case for leptogenesis. 
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or even falsify it, mainly by establishing that the seesaw mechanism is not responsible for the observed 
neutrino masses. By itself, failure in revealing signals of Ouf3f3 decays will not disprove leptogenesis. 
Indeed, with normal neutrino mass hierarchy one expects that the rates of lepton-number-violating 
processes are below experimental sensitivity. However, if neutrinos masses are quasi-degenerate or 
inversely hierarchical, and future measurements of the oscillation parameters will not fluctuate too 
much away from the present best fit values, the most sensitive decay experiments scheduled for 

the near future should be able to detect a signal |251j . If instead the limit on is pushed below 

~ lOmeV (a quite challenging task), this would suggest that either the mass hierarchy is normal, 
or neutrinos are not Majorana particles. The latter possibility would disprove the seesaw model 
and standard leptogenesis. Thus, determining the order of the neutrino mass spectrum is extremely 
important to shed light on the connection between Oi^/3/3 decay experiments and leptogenesis. In 
summary, if it is established that the neutrino mass hierarchy is inverted and at the same time no 
signal of 0iyf3f3 decays is detected at a level |mee| ^ lOmeV, one could conclude that the seesaw is not 
at the origin of the neutrino masses, and that (standard) leptogenesis is not the correct explanation 
of the baryon asymmetry. As concerns CP violation, a failure in detecting leptonic CP violation will 
not weaken the case for leptogenesis in a significant way. Instead, it would mean that the Dirac CP 
phase is small enough to render L-conserving CP-violating effects unobservable. 

Finally, the CERN LHC has the capability of providing information that is relevant to leptogenesis, 
since it can play a fundamental role in establishing that the origin of the neutrino masses is not due 
to the seesaw mechanism, thus leaving no strong motivation for leptogenesis. This may happen in 
several different ways. For example (assuming that the related new physics is discovered), the LHC 
will be able to test if the detailed phenomenology of any of the following models is compatible with 
an explanation of the observed pattern of neutrino masses and mixing angles: supersymmetric R- 
parity violating couplings and/or L- violating bilinear terms \252\ 1253] : leptoquarks [254', '255J; triplet 
Higgses \256i i257j : new scalar particles of the type predicted in the Zee-Babu [ 258, .259j types of 
models \260\ 12611 I262j . It is conceivable that such discoveries can eventually exclude the seesaw 
mechanism and rule out leptogenesis. 
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